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In  this  paper,  we  discuss  the  detection  of  systematic  biases  in  star  positions  of  the  USNO  A1.0,  A2.0,  and 
B1.0  catalogs,  as  deduced  from  the  residuals  of  numbered  asteroid  observations.  We  present  a  technique 
for  the  removal  of  these  biases,  and  validate  this  technique  by  illustrating  the  resulting  improvements  in 
numbered  asteroid  residuals,  and  by  establishing  that  debiased  orbits  predict  omitted  observations  more 
accurately  than  do  orbits  derived  from  non-debiased  observations.  We  also  illustrate  the  benefits  of  debi¬ 
asing  to  high-precision  astrometric  applications  such  as  asteroid  mass  determination  and  collision  anal¬ 
ysis,  including  a  refined  prediction  of  the  impact  probability  of  99942  Apophis.  Specifically,  we  find  the  IP 
of  Apophis  to  be  lowered  by  nearly  an  order  of  magnitude  to  4.5  x  10”6  for  the  2036  close  approach. 
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1.  Introduction 

While  future  asteroid  surveys,  such  as  PanSTARRS,  promise  to 
routinely  provide  observations  with  uncertainties  on  the  order  of 
0.1  arcsec  (Jedicke  et  al.,  2004),  many  high-precision  astrometric 
applications,  such  as  mass  determination  and  Yarkovsky  analysis, 
study  interactions  fixed  in  the  past,  or  rely  upon  long  observational 
baselines  to  model  slowly  evolving  phenomena.  So  while  new 
observations  may  be  helpful,  we  must  also  use  contemporaneous 
and  historical  observations.  Our  long-term  goal,  therefore,  is  to  cre¬ 
ate  a  statistical  error  model  of  asteroid  observations,  providing 
realistic,  observatory-specific  estimates  of  error  correlations  and 
uncertainties  that  will  allow  us  to  make  the  best  possible  use  of 
the  existing  body  of  observations. 

The  astrometric  reduction  of  such  observations  relies  heavily  on 
accurate  star  catalogs;  background  stars  in  an  image  are  identified, 
and  used  as  references  against  which  the  positions  of  the  head  and 
tail  of  an  asteroid  trail  may  be  determined.  Ideally,  star  catalogs 
should  be  “dense  and  deep”,  containing  a  great  many  stars  of  vary¬ 
ing  brightness  distributed  throughout  the  entire  sky,  allowing 
high-precision  astrometry  regardless  of  asteroid  size  or  location. 

In  that  context,  three  all-sky  US  Naval  Observatory  star  catalogs 
are  among  the  most  useful.  The  USNO  A1.0  catalog,  introduced  in 
1996,  contains  488,006,860  sources  down  to  V  magnitude  20,  with 
an  estimated  accuracy  of  0.25  arcsec  (Monet  et  al.,  1998).  The  A2.0 
catalog,  introduced  in  1998,  is  an  update  of  A1.0  that  moves  the 


*  Corresponding  author. 

E-mail  addresses:  steve.chesley@jpl.nasa.gov  (S.R.  Chesley),  jimbaerl@earthlink. 
net  (J.  Baer). 

0019-1035/$  -  see  front  matter  ©  2010  Elsevier  Inc.  All  rights  reserved. 
doi:l  0.1 016/j.icarus.201 0.06.003 


reference  frame  from  the  Guide  Star  Catalog  to  the  ICRF;  it  contains 
526,280,881  sources  (Monet,  1998).  The  B1.0  catalog,  introduced  in 
2003,  contains  1,042,618,261  objects  down  to  V  magnitude  21, 
with  an  estimated  accuracy  of  0.20  arcsec  (Monet  et  al.,  2003). 

Note  that  other  catalogs,  such  as  UCAC2  and  2MASS,  have  smal¬ 
ler  estimated  position  errors;  but  these  catalogs  each  have  limita¬ 
tions  relative  to  B1.0.  The  UCAC2  catalog,  introduced  in  2004, 
contains  48,330,571  stars  with  positions  accurate  to  within 
0.07  arcsec  at  the  limiting  magnitude  of  16;  but  it  only  covers 
the  sky  from  declination  -90°  to  +40°  (Zacharias  et  al.,  2004). 
(Note;  The  recently  introduced  UCAC3  catalog  (Zacharias  et  al., 
2004)  covers  the  entire  sky,  thus  resolving  this  limitation.)  The 
all-sky  2MASS  catalog,  introduced  in  2003,  contains  470,992,970 
objects,  with  positions  accurate  to  within  0.07  arcsec  (Skrutskie 
et  al.,  2006);  but  it  is  limited  to  V  magnitude  17,  meaning  that  ref¬ 
erence  stars  may  be  overexposed  in  images  of  small  asteroids. 

While  creating  the  first  iteration  of  a  statistical  error  model  of 
asteroid  observations,  we  detected  significant  biases  in  the  residu¬ 
als  of  the  numbered  asteroids,  which  were  traced  to  biases  in  the 
star  catalogs  from  which  those  observations  were  reduced. 

After  describing  our  discovery  and  investigation  of  these  biases 
(Section  2),  we  will  describe  the  available  asteroid  astrometric  data 
(Section  3).  Next,  we  present  a  technique  to  debias  these  observa¬ 
tions  (Section  4),  and  demonstrate  that  it  substantially  eliminates 
systematic  errors  in  orbit  fits  of  the  numbered  asteroids  (Section  5). 
We  validate  this  new  debiasing  technique  by  demonstrating  that  it 
produces  orbits  that  predict  better,  which  is  the  crucial  test  of  the 
quality  of  an  orbit  estimate  (Section  6).  Finally,  we  illustrate  the 
application  of  catalog  debiasing  to  issues  such  as  estimating  the 
mass  of  perturbing  asteroids,  and  assessing  the  probability  of  an 
impact  for  potentially  hazardous  asteroids  (Section  7). 
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With  the  catalog  biases  resolved,  a  future  paper  will  describe  the 
subsequent  development  of  a  statistical  error  model  for  asteroid 
observations,  demonstrate  its  validity,  and  detail  its  applications. 

2.  Evidence  for  bias 

Star  catalog  biases  were  encountered  during  the  development 
of  a  statistical  error  model  for  astrometric  asteroid  observations, 
as  described  in  Carpino  et  al.  (2003).  Such  a  model  requires  that 
the  observational  residuals  be  unbiased,  normally  distributed, 
and  uncorrelated. 

In  testing  the  first  iteration  of  our  error  model,  however,  it  be¬ 
came  clear  that  these  assumptions  were  being  violated. 

As  illustrated  in  Fig.  1,  for  instance,  probability  distributions  of 
the  nominal  declination  residuals  for  the  numbered  asteroids  re¬ 
vealed  a  clear  bias  of  approximately  +0.16  arcsec.  Significantly, 
no  such  bias  is  evident  in  the  right  ascension  residuals;  this  led 
us  initially  to  believe  that  only  the  declination  observations  were 
affected. 

Additionally,  we  calculated  a  kurtosis  of  approximately  4.6  for 
both  the  RA  and  the  DEC  residual  probability  distributions  depicted 
in  Fig.  1 ;  this  compares  to  a  kurtosis  of  3.0  for  the  standard  normal 
distribution,  indicating  that  the  RA  and  DEC  distributions  are 
slightly  peaked,  thus  deviating  somewhat  from  the  expected  nor¬ 
mal  distribution.  Adjusting  the  chi-square  observation  rejection 
threshold  in  the  orbit  determination  algorithm  failed  to  reduce 
either  kurtosis  to  expected  levels.  (Note:  The  definition  of  kurtosis 
used  in  this  paper  results  in  a  normal  distribution  having  a  kurtosis 
of  3.0.) 

Finally,  as  will  be  discussed  in  Section  5.3,  the  residuals  from 
closely-spaced  observations  of  the  same  asteroid  made  by  the 
same  observatory  appeared  highly-correlated.  While  the  correla¬ 
tions  associated  with  each  specific  observatory  differed  somewhat 
in  magnitude,  they  remained  significant  even  for  observations  sep¬ 
arated  by  several  days. 

Based  on  these  data,  we  hypothesized  that  there  were  system¬ 
atic  declination  biases  in  the  star  catalogs  used  to  reduce  the  aster¬ 
oid  observations.  And  as  we  investigated  further,  we  encountered 
evidence  outside  of  our  own  work. 

Fig.  2,  dating  to  mid-2008,  depicts  histograms  of  the  postfit 
means  of  right  ascension  and  declination  residuals  for  the  1649 
numbered  asteroids  under  automated  orbit  maintenance  at  that 
time.  While  the  right  ascension  histogram  shows  a  symmetric  dis¬ 
tribution  with  a  mean  of  0.004  arcsec  and  a  standard  deviation  of 


Raw  RA  and  Dec  Residuals 


Fig.  1.  Nominal  residuals  of  the  numbered  asteroids,  illustrating  DEC  bias. 


Mean  of  Normalized  RA  Residuals 


Fig.  2.  Mean  of  normalized  postfit  residuals  for  1649  numbered  asteroids  under 
automated  orbit  maintenance  as  of  mid-2008,  including  all  numbered  NEAs,  plus  a 
handful  of  other  targets  of  interest,  such  as  space  mission  targets.  The  residuals 
were  normalized  by  dividing  each  residual  by  its  corresponding  weight;  and  since 
the  vast  majority  of  observations  are  weighted  at  1  arcsec,  the  abscissa  is 
approximately  in  arcsec. 

0.057  arcsec,  the  declination  histogram  is  decidedly  nonGaussian 
and  asymmetric,  with  a  mean  of  0.074  arcsec  and  a  standard  devi¬ 
ation  of  0.095  arcsec. 

Upon  analysis  of  the  Rosetta  spacecraft’s  encounter  with  Aster¬ 
oid  2867  Steins,  Morley  (T.  Morley,  private  communication)  inde¬ 
pendently  deduced  a  bias  of  +0.212  arcsec  in  the  declination 
observations  of  that  asteroid.  However,  to  fully  account  for  the  ob¬ 
served  encounter  geometry,  a  right  ascension  bias  of  +0.092  arcsec 
was  also  necessary.  This  was  the  first  indication  that  a  bias  was 
present  in  both  coordinates.  As  Morley  noted,  a  systematic  bias 
in  the  right  ascension  observations  of  asteroids  would  not  be 
immediately  obvious,  since  a  least-squares  orbit  determination 
algorithm  would  simply  rotate  the  calculated  orbit  in  the  equato¬ 
rial  frame  about  the  celestial  pole  to  minimize  the  RMS  error,  thus 
eliminating  the  evidence  of  a  net  right  ascension  bias.  Since  none  of 
the  other  orbital  parameters  could  be  manipulated  so  as  to  elimi¬ 
nate  the  declination  bias,  it  appeared  to  be  the  only  systematic  er¬ 
ror  in  the  observations. 

As  will  be  described  more  fully  in  Section  7.2.1,  Mauna  Kea 
Observatory  (MPC  observatory  code  568)  independently  noted  a 
persistent  positive  declination  bias  in  observations  of  99942  Apo- 
phis  (D.  Tholen,  private  communication).  In  the  discussions  that 
followed,  Tholen  referred  us  to  da  Silva  Neto  et  al.  (2005),  who 
had  observed  a  mean  declination  bias  of  approximately  0.11  arcsec 
in  the  USNO  B1.0  positions  of  the  optical  counterparts  of  64  1CRF 
sources.  Using  the  more  accurate  UCAC2  catalog  as  a  reference, 
comparisons  of  the  positions  of  stars  appearing  in  both  B1.0  and 
UCAC2  were  used  to  derive  local  corrections  to  the  B1.0  positions; 
implementing  these  local  corrections,  the  B1.0  declination  biases 
for  the  ICRF  sources  were  reduced  to  approximately  0.03  arcsec. 
In  time,  this  idea  proved  most  useful. 

3.  The  asteroid  astrometric  data 

The  Minor  Planet  Center  (MPC),  hosted  by  the  Smithsonian 
Astrophysical  Observatory  at  Harvard  Univ.,  Cambridge,  Mass., 
operates  under  the  auspices  of  Division  III  of  the  International 
Astronomical  Union  as  the  clearinghouse  for  asteroid  astrometric 
measurements,  both  optical  and  radar  (Williams,  2009).  As  a  part 
of  this  function,  the  MPC  archives  and  distributes  the  available  data 
on  intervals  of  approximately  one  month. 
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Table  1 

Optical  observation  data  types. 


Type 

MPC  flag 

Count 

CCD 

c,  c 

58,083,989 

Former  B1 950.0 

A 

662,222 

Photographic 

303,325 

Transit  circle 

T 

26,968 

Micrometer 

M 

12,081 

Approximate 

X 

10,763 

Hipparcos 

H 

5494 

Occultation 

E 

991 

Satellite 

S 

277 

Roving  observer 

V 

177 

Mis-tagged  obs. 

b 

1 

Total 

- 

59,106,288 

We  have  based  the  present  study  on  the  complete  MPC  astro- 
metric  data  set  dated  2008  December  12.  In  this  data  set  there 
are  59,106,288  optical  observations  of  asteroids,  a  number  that  in¬ 
cludes  48,807,317  optical  observations  of  202,855  numbered  aster¬ 
oids  and  an  additional  10,298,971  optical  observations  of  asteroids 
temporarily  designated  but  not  yet  numbered.  There  are  also  radar 
delay  and  Doppler  observations,  which  we  neglect  for  the  present 
discussion.  The  optical  data  set  is  dominated  by  CCD  observations, 
which  comprise  98.3%  of  the  observations. 

The  MPC  astrometric  data  are  stored  in  80-column  ASCII  re¬ 
cords.  A  nearly  complete  description  of  the  data  format  is  available 
on  the  MPC  website,1  which  we  briefly  summarize  here.  There  is  one 
record  per  observation,  except  for  Satellite  and  Roving  observations, 
which  require  a  second  observation  record  to  describe  the  position  of 
the  observer.  The  primary  information  on  the  optical  astrometry 
data  record  is  the  observation  time  and  sky  position,  given  as 
(J2000)  right  ascension  and  declination.  The  measurement  tech¬ 
nique,  or  observation  type,  is  recorded  by  a  code  in  column  15,  the 
counts  for  which  are  given  in  Table  1.  The  observatory  code  in  col¬ 
umns  78-80  is  a  three-character  identifier  that  uniquely  identifies 
the  observing  location,  and  sometimes  the  observing  program.  (In 
cases  where  multiple  observing  programs  share  the  same  observa¬ 
tory  code,  e.g.,  Mauna  Kea  or  Mt.  Palomar,  there  may  be  an  addi¬ 
tional  field  that  records  the  observing  program.)  The  observatory 
codes  and  the  names  and  locations  of  the  associated  observatories 
are  tabulated  by  the  MPC.2  In  this  paper,  for  brevity,  we  will  often 
refer  to  the  observatories  only  by  their  observatory  code.  There  is 
also  a  field  for  the  observation’s  Minor  Planet  Circular  or  Supplement 
reference,  and  a  field  that  allows  the  observer  to  provide  a  terse 
annotation  regarding  the  observation  quality  or  reduction  technique. 
Photometry  information  may  also  be  part  of  the  record. 

For  the  present  effort,  the  most  important  part  of  the  observa¬ 
tion  record  is  the  undocumented  annotation,  in  column  72,  of 
the  star  catalog  that  was  used  in  the  reduction  of  the  astrometry 
(Gareth  Williams,  private  communication).  Table  2  lists  the  field 
values  and  their  meanings,  along  with  the  number  of  CCD  observa¬ 
tions  for  each  field  value.  An  inspection  of  the  data  indicates  that 
the  MPC  did  not  start  recording  the  catalog  in  the  observation  re¬ 
cord  until  mid-2001,  and  even  after  that  many  observers  did  not 
report  the  catalog  information  to  the  MPC.  Thus  for  about  17%  of 
the  astrometry  there  is  no  catalog  information  available  on  the 
MPC  data  record.  Fortunately,  much  of  the  missing  catalog  infor¬ 
mation  is  from  only  a  few  observatory  codes,  and  for  data  from 
prior  to  mid-2001.  Thus,  although  the  MPC  did  not  record  the 
information,  some  of  the  observing  programs  still  have  their  obser¬ 
vation  logs  from  that  period  and  so  we  were  able  to  supply  catalog 


1  http://www.cfa.harvard.edu/iau/info/ObsFormat.html. 

2  http://www.cfa.harvard.edu/iau/lists/ObsCodesF.html. 


Table  2 

Catalog  flags  found  for  CCD  observations  on  MPC  data 
file. 


Catalog 

MPC  flag 

Count 

USNO  A1.0 

a 

173,556 

USNO  SA1.0 

b 

31,044 

USNO  A2.0 

c 

26,071,594 

USNO  SA2.0 

d 

1,327,783 

UCAC1 

e 

243,111 

Tycho- 1 

f 

0 

Tycho-2 

g 

299,648 

GSC  1.0 

h 

0 

GSC  1.1 

i 

10,055 

GSC  1.2 

j 

11,665 

GSC  2.2 

k 

180 

ACT 

1 

13,866 

GSC  ACT 

m 

364,575 

TRC 

n 

0 

USNO  B1.0 

0 

6,799,942 

PPM 

p 

0 

UCAC2-beta 

q 

0 

UCAC2 

r 

12,315,889 

USNO  B2.0 

s 

445 

UCAC3-beta 

t 

0 

UCAC3 

u 

0 

NOMAD 

V 

2763 

CMC 

w 

337,053 

Hip  2 

X 

0 

GSC  (generic) 

z 

1537 

Unspecified 

- 

10,079,283 

Total 

- 

58,083,989 

Table  3 

Observer-supplied  catalog  information. 


Obs.  code 

Time  frame 

Catalog 

704 

Before  2000-January-01 

USNO  A1.0 

704 

After  2000-January-01 

USNO  A2.0 

691 

1 991  -August-3 1  -1 999-September-28 

GSC  1 

691 

1999-September-29-2000-December-21 

USNO  A1.0 

691 

2000-December-22-2006-December-26 

USNO  A2.0 

703 

Before  2005-January-01 

USNO  A2.0 

info  for  59%  of  the  untagged  astrometry.  Table  3  lists  the  rules  that 
we  have  applied  to  otherwise  untagged  astrometry. 

Many  of  the  individual  catalog  tags  given  in  Table  2  can  be 
grouped  together  for  our  purposes.  For  instance,  the  USNO  SA2.0 
catalog  is  merely  a  subset  of  USNO  A2.0,  and  UCAC2  is  an  extension 
to  UCAC1.  We  have  merged  the  various  column-72  symbols  as  de¬ 
tailed  by  Table  4.  We  note  that  over  90%  of  the  astrometry  is  com¬ 
prised  from  four  catalogs,  USNO  A1.0,  USNO  A2.0,  USNO  B1.0,  and 
UCAC.  An  additional  7.2%  still  have  catalog  unknown,  although 
most  of  these  are  from  major  NEO  survey  programs,  and  so  these 
also  are  probably  dominated  by  the  same  four  catalogs.  As  we  dis¬ 
cuss  below,  one  can  infer  which  catalog  was  used  in  such  cases  by  a 
statistical  analysis  of  the  residuals.  Fig.  3  depicts  the  history  of  cat¬ 
alog  use  for  the  data  given  in  Table  4.  Some  catalogs  appear  earlier 
than  their  release  date  because  older  images  were  in  some  cases 
remeasured  with  new  software  and  star  catalogs.  USNO  A2.0  has 
been  used  heavily,  with  more  than  2  million  observations  every 
year  since  2000,  which  happened  to  be  the  last  year  that  USNO 
A1.0  saw  extensive  usage.  But  starting  in  2005,  the  UCAC  and 
USNO  B1 .0  catalogs  have  together  accounted  for  the  majority  of  re¬ 
ported  astrometry. 

4.  Two  approaches  to  debiasing 

There  are  two  obvious  way  to  quantify  the  effect  of  star  catalog 
bias  in  asteroid  orbit  determination.  One  can  proceed  indirectly, 
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Table  4 

Merged  catalog  counts. 


Catalog 

MPC  flags 

Count 

% 

USNO  A2.0 

c,  d 

30,786,427 

53.0 

UCAC 

e,  q,  r,  t,  u 

12,559,000 

21.6 

USNO  B1.0 

o,  s 

6,800,387 

11.7 

Unknown 

- 

4,171,601 

7.2 

USNO  A1.0 

a,  b 

2,205,452 

3.8 

CSC  1 

h,  i,  j,  z 

543,037 

0.9 

CSC  ACT 

m 

364,575 

0.6 

CMC 

w 

337,053 

0.6 

Tycho 

f.  g 

299,648 

0.5 

ACT 

i 

13,866 

0.0 

NOMAD 

V 

2,763 

0.0 

CSC  2 

k 

180 

0.0 

Total 

- 

58,083,989 

- 

through  a  systematic  study  of  the  postfit  residuals  of  astrometric 
asteroid  observations,  or  directly,  by  comparing  star  catalog  posi¬ 
tions  with  those  in  a  more  accurate  and  precise  reference.  Since 
the  problem  first  arose  for  us  from  the  perspective  of  asteroid 
residuals,  we  began  with  the  first  approach. 

4.1.  Debiasing  from  indirect  reference  to  postfit  astrometric  residuals 

As  shown  in  Fig.  1,  the  observations  of  numbered  asteroids 
demonstrated  a  clear  observational  bias  in  declination;  but  no  such 
bias  appeared  in  right  ascension.  We  therefore  began  by  debiasing 
only  the  declination  coordinates.  After  calculating  the  orbits  of  the 
numbered  asteroids  and  compiling  the  “observed-predicted” 
residuals,  we  divided  the  sky  into  18  zones  of  declination,  each 
10°  wide;  within  each  zone,  we  calculated  the  mean  declination 
biases  for  USNO  A1.0,  A2.0,  B1.0,  UCAC2,  and  observations  of  vari¬ 
ous  eras  whose  reduction  catalogs  were  unknown.  We  then  cre¬ 
ated  best-fit  polynomials  for  each  of  these  groups  to  model  the 
declination  biases  in  each  zone,  and  used  them  to  debias  the  obser¬ 
vations  for  the  next  iteration. 

Initially,  the  approach  seemed  promising;  the  biases  resulting 
from  the  second  iteration  were  reduced  to  one  third  of  their  previ¬ 
ous  values.  However,  in  order  to  achieve  these  reductions,  we  had 
to  apply  significant  corrections  (on  the  order  of  0.35  arcsec)  to  the 
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Fig.  3.  Year-by-year  counts  of  star  catalog  usage  for  CCD  astrometry. 


UCAC2  declination  observations.  Given  the  stated  accuracy  of  the 
UCAC2  star  positions,  this  was  clearly  unrealistic.  We  concluded 
that  the  biases  in  the  other  catalogs  were  likely  resulting  in  best- 
fit  orbits  that  made  the  UCAC2  residuals  appear  much  larger  than 
their  intrinsic  errors.  Therefore,  we  decided  to  repeat  the  second 
iteration,  but  with  no  corrections  applied  to  the  UCAC2  derived 
positions. 

The  results  from  the  revised  second  iteration  showed  significant 
improvement;  for  most  of  the  observation  sets,  the  biases  had  been 
reduced  to  below  0.06  arcsec.  And  yet  the  biases  in  the  UCAC2  de¬ 
rived  observations  remained  above  0.2  arcsec.  A  third  iteration 
yielded  little  further  reduction. 

By  this  point,  additional  analysis  had  demonstrated  that  the 
asteroid  residuals  did  indeed  contain  a  right  ascension  bias;  and 
we  suspected  that  our  failure  to  account  for  this  was  responsible 
for  the  slow  convergence  in  eliminating  the  observed  bias  in  the 
UCAC2  residuals.  We  therefore  implemented  a  two-dimensional 
bias  model  for  all  future  iterations,  with  a  grid  spacing  of  1  h  in 
right  ascension  and  10°  in  declination.  As  before,  the  mean  ob¬ 
served  bias  in  each  grid  square  was  used  to  debias  the  observations 
in  the  next  iteration.  But  unfortunately,  two  further  iterations 
yielded  little  further  improvement;  the  biases  in  the  UCAC2  de¬ 
rived  observations  remained  near  0.2  arcsec. 

Despite  these  difficulties,  the  bias  corrections  that  we  derived 
through  this  indirect  approach  do  share  many  of  the  gross  proper¬ 
ties  seen  in  the  corresponding  (and  much  more  detailed)  sky  maps 
derived  from  the  direct  approach  seen  in  the  next  section.  In  retro¬ 
spect,  we  believe  that,  had  we  started  with  a  two-dimensional  bias 
model  in  the  first  place,  perhaps  with  smaller  bins  on  the  sky,  this 
iterative  technique  might  have  yielded  adequate  results;  but  we  do 
not  believe  that  it  can  be  as  effective  as  the  direct  approach  that  we 
have  developed,  which  we  describe  now. 

4.2.  Debiasing  from  direct  reference  to  star  catalogs 

An  assessment  of  the  systematic  errors  in  star  catalog  astrome¬ 
try  is  the  key  ingredient  in  the  direct  approach  to  quantification  of 
systematic  errors  in  asteroid  astrometry  from  this  source.  Thus  we 
needed  to  produce  a  spatially  resolved  astrometric  comparison  of 
various  catalogs  that  have  been  used  as  reference  catalogs  for  the 
astrometry  of  solar  system  objects.  Because  of  its  wide  dynamic 
range  and  its  rigorous  ties  to  the  ICRS  (Arias  et  al.,  1995)  using 
the  Tycho-2  catalog  (Hog  et  al.,  2000),  the  2MASS  star  catalog  (Skr- 
utskie  et  al„  2006)  was  chosen  as  the  catalog  to  which  all  others 
were  compared.  According  to  its  documentation,  the  2MASS  sys¬ 
tematic  position  errors  are  70-80  milliarcsec  with  respect  to  the 
ICRS.  The  2MASS  catalog  does  not  contain  proper  motions,  but 
the  mean  epoch  of  the  2MASS  observations  is  very  close  to 
2000.0.  Hence,  the  comparison  of  all  of  the  catalogs  ignored  proper 
motions. 

To  achieve  the  desired  spatial  resolution,  the  JPL  HEALPix3  tes¬ 
sellation  (Gorski  et  al.,  2005)  was  used.  Adopting  SIDE  =  64  with 
the  HEALPix  package  produces  49,152  tiles  on  the  sky,  called  HEAL¬ 
Pix,  each  with  an  area  just  a  bit  smaller  than  one  square  degree.  This 
seems  to  be  a  reasonable  compromise  between  high  resolution  and 
having  enough  stars  in  each  tile.  Given  this  tessellation  scheme,  the 
tile  number  for  each  entry  of  each  star  catalog  was  computed,  and 
the  catalogs  were  sorted  on  this  index  to  expedite  the  calculation. 

The  numerical  processing  was  the  same  for  all  catalogs  under 
consideration.  For  each  HEALPix,  the  stars  in  the  catalog  and  in 
2MASS  were  extracted,  and  a  spatial  correlation  using  a  radius  of 
2.0  arcsec  was  used  to  identify  stars  in  common.  The  size  of  the  ra¬ 
dius  was  chosen  to  be  large  enough  to  accommodate  the  expected 
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systematic  differences  (up  to  1.0  arcsec),  but  small  enough  to  min¬ 
imize  the  number  of  spurious  correlations.  Once  the  pairs  were 
identified,  the  mean  and  standard  deviation  in  the  right  ascension 
and  declination  directions  were  computed.  These  provide  a  mea¬ 
sure  of  the  systematic  and  random  local  position  errors.  Note  that 
the  calculation  was  done  in  the  tangent  plane.  For  each  HEALPix,  a 
tangent  point  at  the  center  of  the  tile  was  adopted,  and  the  tangent 
plane  positions  for  the  catalog  and  2MASS  stars  were  computed. 
Hence,  the  resulting  astrometric  shifts  and  dispersion  should  mi¬ 
mic  those  needed  to  correct  the  analysis  of  observational  data. 

We  performed  the  analysis  for  five  of  the  catalogs  listed  in 
Table  4:  USNO  A1.0,  USNO  A2.0,  USNO  B1.0,  UCAC2,  and  Tycho. 
The  output  of  the  comparison  of  these  catalogs  with  2MASS  was 
a  single  file  containing  49,152  lines  each  with  the  nominal  position 
of  the  HEALPix  on  the  sky  and  values  for  the  mean  and  standard 
deviation  for  the  astrometric  offset  (excepting  for  UCAC2  which 
does  not  cover  the  entire  sky).  While  the  primary  use  of  the  com¬ 
parison  file  is  intended  to  be  a  scheme  by  which  individual  obser¬ 
vations  can  be  corrected  to  the  astrometric  reference  frame  defined 
by  the  2MASS  catalog,  various  visualization  tools  were  used  to 
understand  the  gross  properties  of  the  various  catalogs. 

Fig.  4  includes  plots  of  the  probability  density  of  the  mean  dif¬ 
ferences  between  the  various  catalogs  and  2MASS  among  the 
49,1 52  HEALPix  cells.  Table  5  lists  the  mean  and  standard  deviation 
of  the  same  distributions.  From  these,  we  observe  that  the  preci¬ 
sion  and  accuracy  of  the  Tycho  and  UCAC  catalogs  are  comparable, 
and  vastly  superior  to  the  other  catalogs.  Indeed,  the  differences 
with  respect  to  2MASS  appear  small  enough  that  it  is  difficult  to 
know  whether  they  should  be  attributed  to  the  subject  catalogs, 
to  2MASS  or  whether  the  differences  are  a  blend  of  errors  from 
both  sides  of  the  comparison. 

But,  Tycho  and  UCAC  have  relatively  few  stars  and  are  too 
sparse  for  most  asteroid  astrometry  work,  and  so  much  of  the 
asteroid  astrometry  has  been  and  continues  to  be  measured 
against  the  much  more  densely  populated  USNO  family  of  catalogs. 
Here  we  see  that  the  USNO  A1.0  has  relatively  large  dispersions,  of 
order  400  mas,  but  overall  it  exhibits  a  relatively  modest  bias.  On 
the  other  hand,  the  USNO  A2.0  and  B1.0  catalogs  reveal  significant 
biases  overall,  especially  in  declination,  with  mean  offsets  of  142 
mas  and  126  mas,  respectively.  These  values  are  comparable  to 
the  standard  deviations  of  the  samples.  In  right  ascension,  B1.0  ap¬ 
pears  relatively  unbiased  overall,  while  A2.0  does  show  a  modest 
net  RA  bias. 
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Fig.  4.  Probability  densities  of  inter-catalog  systematic  errors,  as  compared  to  the 
2MASS  catalog.  For  each  plot,  the  abscissa  is  the  difference  in  arcsec  between  the 
given  catalog  and  2MASS  and  the  ordinate  is  the  associated  probability  density  in 
arcsec-1.  Note  that  the  plots  are  not  all  on  the  same  scale. 


Table  5 

Inter-catalog  systematic  errors,  with  respect  to  2MASS. 


Catalog 

Mean  ±  std.  dev. 

RA  (mas) 

DEC  (mas) 

Tycho-2 

-1  ±28 

-12  ±  24 

UCAC2 

2  ±23 

-7  ±20 

USNO  B1.0 

-16  ±  123 

126  ±123 

USNO  A2.0 

63  ±180 

142  ±189 

USNO  A1.0 

-41  ±419 

-34  ±352 

Of  course,  the  discussion  so  far  does  not  reveal  anything  about 
the  structure  on  the  sky.  If  half  of  the  sky  has  a  positive  offset  and 
half  of  the  sky  has  a  negative  offset  then  the  mean  offset  could  still 
be  small.  To  address  this  question,  Fig.  5  depicts  sky  maps  of  the  RA 
and  DEC  differences  with  respect  to  2MASS  for  the  five  catalogs  un¬ 
der  consideration;  while  Fig.  6  depicts  the  standard  deviations  in 
those  RA  differences.  Note  first  that  the  color  scale  is  different  for 
different  catalogs  in  order  to  show  as  much  of  the  structure  in 
the  catalog  bias  as  possible.  In  these  plots,  regions  with  relatively 
small  biases  are  colored  green,  while  red  regions  show  positive 
bias  and  blue  regions  have  negative  bias.  The  USNO  A1.0  catalog 
does  show  significant  structure,  despite  its  rather  innocuous  pre¬ 
sentation  in  Fig.  4,  and  yet  the  hemispheric  differences  are  rela¬ 
tively  modest.  On  the  other  hand,  the  USNO  A2.0  and  B1.0 
catalogs,  which  together  account  for  approximately  70%  of  all 
asteroid  astrometry,  are  more  problematic.  We  note  in  particular 
that  the  declination  maps  for  these  two  catalogs  are  plainly  “hot” 
overall,  in  agreement  with  Fig.  4.  Also,  the  RA  maps  show  evidence 
for  hemispheric  bias,  with  a  large  cold  region  (negative  offset)  cen¬ 
tered  around  45°  RA  extending  from  pole  to  pole,  and  an  even 
stronger  hot  region  opposite,  which  is,  however,  roughly  limited 
to  the  northern  hemisphere  for  A2.0  and  the  southern  hemisphere 
for  B1 .0.  In  contrast,  the  UCAC  and  Tycho  sky  maps  show  very  lim¬ 
ited  structure,  more  consistent  with  random  fluctuations. 

The  procedure  followed  in  this  paper  is  to  use  the  mean  catalog 
differences  on  each  HEALPix  cell  to  debias  astrometry  from  the  ass- 
sociated  catalog  in  that  cell.  Specifically,  for  a  given  astrometric 
observation  with  known  catalog,  we  first  compute  which  of  the 
49,152  HEALPix  sky  tiles  contains  the  astrometry  and  look  up  the 
associated  RA  and  DEC  offsets,  which  are  subtracted  from  the 
raw  astrometry  to  produce  debiased  astrometry.  In  light  of  the  fore¬ 
going  discussion  on  the  low  systematic  errors  seen  in  the  UCAC  and 
Tycho  catalogs,  we  consider  that  astrometry  reduced  with  these 
catalogs  should  not  show  significant  catalog-based  systematic  bias. 
Therefore,  we  apply  catalog  bias  corrections  only  for  astrometry 
derived  from  the  USNO  A1.0,  A2.0  and  B1.0  catalogs.  The  catalog 
comparison  file  and  ancillary  documentation  are  available  to  other 
researchers  via  FTP.4 

5.  A  test  on  all  numbered  asteroids 

To  evaluate  the  effectiveness  of  the  debiasing  technique  out¬ 
lined  in  the  previous  section,  we  have  calculated  the  orbits  of 
asteroids,  both  without  and  with  debiasing  corrections,  to  obtain 
what  we  call  "raw"  and  “debiased"  fits,  respectively.  We  limit  the 
objects  considered  to  those  with  well-determined  orbits,  namely 
the  numbered  asteroids,  to  ensure  that  random  errors  do  not— in 
the  mean— contaminate  the  results  in  a  significant  way.  Also,  be¬ 
cause  the  objective  here  is  not  necessarily  to  obtain  the  best  possi¬ 
ble  orbits,  but  rather  to  evaluate  the  efficacy  of  debiasing,  we  have 
used  only  data  that  comes  within  the  scope  of  the  debiasing  tech¬ 
nique,  namely  astrometry  known  to  be  reduced  with  respect  to  one 
of  the  five  catalogs  discussed  in  the  previous  section.  This  amounts 


4  ftp://sad.jpl.nasa.gov/pub/ssd/debias/debias.tgz. 
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Fig.  5.  HEALPix  maps  of  RA  and  DEC  bias  with  respect  to  2MASS  for  five  star  catalogs,  as  described  in  text. 
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UCAC 


to  48,621,261  observations  of  201,804  asteroids  over  the  interval 
1949-11-19  to  2008-12-12.  (Note:  The  1949-era  observations  are 
likely  precovery  attempts  reduced  using  one  of  the  five  catalogs.) 

Outlier  rejection  was  used  in  the  orbit  calculation,  according  to 
the  technique  described  by  Carpino  et  al.  (2003),  with  separate  fits 
for  the  raw  and  debiased  cases. 

The  distribution  of  the  change  in  orbital  elements  between  the 
raw  and  debiased  fits  is  depicted  in  Fig.  7.  There  we  see  significant 
changes  in  the  inclinations  of  the  debiased  orbits,  relative  to  their 
“raw”  orbit  counterparts.  The  inclinations  generally  declined, 
which  is  consistent  with  the  fact  that  the  majority  of  the  astrome¬ 
try  falls  in  the  north  and  has  a  positive  declination  bias,  which  will 
tend  to  erroneously  increase  the  estimated  inclination.  There  are 
smaller  changes  in  the  semimajor  axis  and  eccentricity  estimates, 
but  still  many  cases  with  >1  sigma  deviations. 

5.1.  Correlation  of  postfit  residuals  with  sky  position 

Fig.  8  depicts  the  postfit  residuals  for  the  raw  fits,  where  the  sky 
map  reflects  the  mean  postfit  residual  in  each  HEALPix  cell.  As  is 
the  case  with  the  asteroids  themselves,  the  majority  of  numbered 
asteroid  detections  are  found  along  the  ecliptic.  Only  cells  with  at 
least  10  residual  measurements  are  plotted.  For  the  USNO  A1.0, 
A2.0  and  B1.0  catalogs,  the  raw  astrometric  residuals  show  a 
remarkably  similar  structure  and  magnitude  to  that  of  the  corre¬ 
sponding  catalog  (Fig.  5).  This  suggests  that,  where  the  observation 
density  is  sufficient,  the  indirect  approach  to  debiasing  (by  refer¬ 
ence  to  postfit  residuals)  may  have  utility  in  cases  where  the  cata¬ 
log  bias  is  poorly  constrained. 

There  is  also  a  significant  bias  in  the  postfit  residuals  of  astrom¬ 
etry  reduced  with  the  Tycho  and  UCAC  catalogs,  despite  the  low 
systematic  errors  seen  with  these  catalogs.  This  is  not  particularly 
surprising,  given  that  the  vast  majority  of  observations  arise  from 
the  significantly  biased  A2.0  and  B1.0  catalogs;  these  catalogs  have 
skewed  the  orbit  estimates  for  the  raw  fits,  leading  to  higher  resid¬ 
uals  in  the  unbiased  Tycho  and  UCAC  astrometry.  In  particular,  po¬ 
sitive  declination  biases  will  tend  to  skew  the  orbital  inclination  so 
that  the  declination  residuals  from  unbiased  observations  will  also 
be  positive. 

The  postfit  residuals  for  the  debiased  orbital  fits  are  depicted  in 
Fig.  9,  where  the  signature  seen  in  Fig.  8  is  now  substantially  re¬ 
duced,  indicating  that  the  biases  have  been  largely  eliminated 
(note,  however,  the  evident  striping  in  the  UCAC  plots,  which  we 
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Fig.  6.  HEALPix  maps  of  l  ight  ascension  standard  deviation  with  respect  to  2MASS.  Fig.  7.  Distribution  of  orbital  element  variation  between  raw  and  debiased  fits  to  all 
Declination  dispersions  are  not  significantly  different.  numbered  asteroids. 
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Fig.  8.  Mean  residuals  for  fits  of  raw  astrometry  to  numbered  asteroids.  The  plots  depict  49,152  equal  area  cells  using  the  HEALPix  algorithm  (Gorski  et  al.,  2005).  Cells  with 
10  or  fewer  observations  are  not  plotted. 
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Fig.  9.  Mean  residuals  for  fits  of  debiased  astrometry  to  numbered  asteroids,  as  in  Fig.  8.  The  debiasing  has  substantially  removed  the  catalog  bias  signal  seen  in  Fig.  8. 
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Fig.  10.  Distribution  of  postfit  residuals  for  raw  and  debiased  orbital  fits  for  all  numbered  asteroids. 


will  discuss  further  in  Section  5.4).  The  histograms  depicted  in 
Fig.  10  indicate  that  the  debiasing  technique  does  indeed  improve 
the  residual  distribution;  in  each  case,  the  curves  representing 
debiased  residuals  have  higher  peaks,  and  are  more  nearly  cen¬ 
tered  on  zero,  indicating  that  the  mean  residuals  and  variances 
have  been  reduced. 

Table  6  quantifies  the  reduction  in  mean  residuals  for  each  cat¬ 
alog,  and  also  indicates  the  scatter  in  the  astrometric  residuals,  be¬ 
fore  and  after  debiasing.  The  table  lists  the  number  of  observations, 
separated  by  catalog  and  observatory,  and  the  associated  RA  and 
DEC  residual  mean  (a ,5)  and  standard  deviation  ((a), (5)).  The  ra¬ 
tios  between  the  debiased  and  raw  numbers  are  also  present. 

Overall,  the  star  catalog-based  debiasing  was  quite  effective, 
with  the  mean  residuals  reduced  by  as  much  as  70%;  but  there 
are  a  few  curious  entries  in  Table  6  that  are  worth  noting.  Of  par¬ 
ticular  interest  are  the  cases  that  still  show  a  moderate  residual 
bias  in  the  debiased  fits.  For  example,  observatory  code  703  (Cata¬ 
lina  Sky  Survey,  Mt.  Bigelow  Schmidt  Camera)  shows  significant  RA 
bias  in  both  the  USNO  A2.0  and  UCAC  astrometry;  the  UCAC  DEC 


bias  is  rather  large,  as  well.  Similarly,  El 2  (Siding  Spring  Survey), 
shows  significant  DEC  bias;  and  a  few  other  observatories  appear 
to  have  some  signature  as  well.  Also,  while  the  star  catalog-based 
debiasing  clearly  reduced  the  mean  residuals  for  observations  re¬ 
duced  using  USNO  B1.0,  the  mean  RA  residuals  for  observatories 
699  and  691  are  still  significant;  in  effect,  their  large  but  opposite 
biases  cancel  each  other  out.  We  propose  a  likely  explanation  for 
these  behaviors  below. 

5.2.  Untagged  astrometry 

Table  6  also  details  the  statistics  of  the  (prefit)  residuals  associ¬ 
ated  with  observations  reduced  using  unknown  catalogs,  which 
appear  dominated  by  observatories  699  (LONEOS)  and  644 
(NEAT-Haleakala).  An  examination  of  the  residuals  allows  us  to  in¬ 
fer  which  catalog  was  used  in  each  case. 

Table  7  contains  a  subset  of  Table  6,  showing  the  mean  residuals 
and  standard  deviations  for  the  five  catalogs,  as  well  as  the  prefit 
residuals  for  the  unknown  observations  from  699,  644,  608,  106, 


Table  6 

Residual  statistics  for  numbered  asteroids  by  catalog  and  observatory. 


Code 

Used 

Percent 

Deleted 

Debiased  (arcsec) 

Raw  (arcsec) 

Ratio  debiased/raw 

a 

(a) 

3 

m 

a 

(a) 

3 

a 

<a) 

3 

w 

USNO  A2.0 

ALL 

26,215,118 

53.93 

211,513 

0.03 

0.61 

0.08 

0.58 

0.06 

0.64 

0.29 

0.60 

0.44 

0.96 

0.27 

0.96 

704 

20,358,631 

77.68 

169,849 

0.05 

0.63 

0.10 

0.60 

0.08 

0.66 

0.30 

0.63 

0.59 

0.96 

0.32 

0.96 

699 

1,158,767 

4.406 

5,577 

0.00 

0.62 

0.11 

0.53 

0.02 

0.65 

0.35 

0.55 

0.08 

0.96 

0.32 

0.95 

691 

1,796,063 

6.812 

4,076 

-0.03 

0.30 

-0.01 

0.30 

0.04 

0.33 

0.23 

0.35 

-0.73 

0.91 

-0.05 

0.84 

608 

576,376 

2.210 

7,582 

0.06 

0.61 

-0.11 

0.75 

0.10 

0.65 

0.13 

0.76 

0.60 

0.95 

-0.86 

0.98 

703 

861,322 

3.317 

15,267 

-0.20 

0.69 

-0.02 

0.63 

-0.15 

0.70 

0.20 

0.65 

1.37 

0.98 

-0.09 

0.97 

644 

707,790 

2.679 

203 

0.00 

0.29 

0.09 

0.30 

0.03 

0.34 

0.29 

0.34 

-0.10 

0.86 

0.29 

0.87 

291 

222,085 

0.842 

359 

-0.10 

0.46 

0.04 

0.32 

-0.03 

0.46 

0.28 

0.36 

2.83 

0.98 

0.14 

0.88 

599 

74,867 

0.285 

469 

-0.01 

0.39 

0.01 

0.34 

0.01 

0.41 

0.26 

0.38 

-1.23 

0.97 

0.05 

0.88 

333 

13,509 

0.051 

76 

0.04 

0.55 

-0.24 

0.53 

0.05 

0.55 

-0.05 

0.54 

0.88 

0.99 

4.63 

0.98 

D35 

29,074 

0.110 

16 

0.11 

0.39 

0.02 

0.38 

0.20 

0.41 

0.21 

0.40 

0.53 

0.95 

0.09 

0.94 

OTH 

416,634 

1.607 

8039 

-0.03 

0.63 

-0.04 

0.51 

0.01 

0.62 

0.16 

0.53 

-1.97 

1.00 

-0.25 

0.96 

UCAC1  &  2 

ALL 

9,109,348 

18.74 

48,420 

-0.10 

0.53 

0.05 

0.49 

-0.13 

0.54 

0.03 

0.52 

0.74 

0.98 

1.91 

0.95 

703 

4,186,043 

46.04 

30,218 

-0.20 

0.63 

0.12 

0.59 

-0.23 

0.64 

0.07 

0.61 

0.84 

0.98 

1.65 

0.98 

G96 

3,283,638 

35.89 

2,974 

-0.01 

0.32 

-0.05 

0.27 

-0.04 

0.33 

-0.09 

0.31 

0.27 

0.95 

0.61 

0.87 

E12 

1,001,493 

10.98 

3,779 

-0.02 

0.50 

0.17 

0.45 

-0.05 

0.50 

0.28 

0.47 

0.46 

0.98 

0.63 

0.96 

683 

200,152 

2.253 

6,218 

0.04 

0.79 

0.07 

0.90 

0.01 

0.80 

0.02 

0.91 

3.32 

0.99 

4.52 

0.99 

J75 

117,112 

1.280 

147 

0.03 

0.41 

-0.06 

0.37 

0.04 

0.42 

-0.10 

0.39 

0.80 

0.98 

0.55 

0.93 

106 

92,049 

1.006 

44 

0.02 

0.40 

-0.05 

0.39 

-0.02 

0.42 

-0.14 

0.41 

-1.01 

0.97 

0.33 

0.95 

143 

1,443 

0.016 

26 

0.09 

0.57 

-0.04 

0.47 

0.07 

0.57 

-0.12 

0.48 

1.26 

1.00 

0.31 

0.97 

OTH 

227,418 

2.538 

5014 

0.01 

0.51 

-0.00 

0.42 

-0.03 

0.51 

-0.04 

0.45 

-0.26 

0.99 

0.02 

0.94 

USNO  B1.0 

ALL 

5,085,749 

10.46 

15,652 

0.03 

0.48 

0.09 

0.42 

-0.06 

0.49 

0.18 

0.44 

-0.40 

0.99 

0.51 

0.95 

699 

2,408,525 

47.44 

11,719 

0.11 

0.61 

0.15 

0.54 

0.02 

0.61 

0.23 

0.56 

6.50 

0.99 

0.63 

0.97 

644 

1,275,002 

25.00 

182 

0.04 

0.24 

0.07 

0.20 

-0.05 

0.26 

0.18 

0.26 

-0.68 

0.92 

0.41 

0.77 

691 

1,166,659 

22.92 

2,366 

-0.16 

0.30 

0.00 

0.28 

-0.23 

0.32 

0.08 

0.32 

0.70 

0.93 

0.06 

0.89 

291 

63,199 

1.239 

12 

-0.22 

0.39 

0.07 

0.26 

-0.29 

0.40 

0.15 

0.29 

0.76 

0.96 

0.49 

0.89 

OTH 

172,364 

3.406 

1373 

0.05 

0.48 

0.03 

0.39 

-0.04 

0.46 

0.14 

0.41 

-1.15 

1.02 

0.24 

0.94 

USNO  A1.0 

ALL 

1,940,649 

3.992 

67,389 

-0.01 

0.72 

0.03 

0.69 

-0.06 

0.76 

-0.03 

0.76 

0.12 

0.95 

-0.93 

0.92 

704 

1,597,775 

82.83 

65,408 

0.00 

0.76 

0.03 

0.73 

-0.05 

0.81 

-0.03 

0.80 

-0.01 

0.95 

-1.06 

0.91 

691 

319,034 

15.97 

1594 

-0.06 

0.49 

0.01 

0.46 

-0.12 

0.51 

-0.06 

0.48 

0.47 

0.96 

-0.22 

0.97 

OTH 

23,840 

1.206 

387 

0.08 

0.65 

0.21 

0.62 

-0.01 

0.63 

0.19 

0.61 

-5.63 

1.04 

1.12 

1.03 

Tycho-2 

ALL 

220,629 

0.46 

564 

0.03 

0.24 

0.01 

0.25 

-0.01 

0.26 

0.00 

0.28 

-3.85 

0.91 

-5.38 

0.89 

689 

133,287 

60.26 

11 

0.00 

0.20 

-0.03 

0.21 

-0.06 

0.22 

-0.04 

0.24 

-0.06 

0.91 

0.83 

0.87 

OTH 

87,342 

39.74 

553 

0.06 

0.29 

0.07 

0.30 

0.07 

0.30 

0.05 

0.34 

0.89 

0.95 

1.29 

0.91 

Unknown 

ALL 

5,336,727 

N/A 

359,619 

0.02 

0.73 

0.16 

0.68 

-0.07 

0.73 

0.13 

0.70 

-0.31 

1.00 

1.23 

0.98 

699 

1,180,015 

20.84 

7,153 

0.11 

0.68 

0.31 

0.57 

0.03 

0.67 

0.27 

0.60 

3.65 

1.02 

1.14 

0.96 

644 

1,119,458 

19.85 

11,277 

0.04 

0.51 

0.43 

0.51 

-0.04 

0.49 

0.40 

0.52 

-1.06 

1.04 

1.08 

0.97 

608 

581,208 

10.43 

13,175 

0.00 

0.78 

0.12 

0.80 

-0.09 

0.77 

0.11 

0.80 

-0.02 

1.02 

1.07 

0.99 

D29 

147,281 

2.591 

325 

0.08 

0.44 

-0.10 

0.41 

0.03 

0.44 

-0.16 

0.44 

2.85 

0.98 

0.60 

0.94 

689 

124,658 

2.194 

330 

-0.02 

0.28 

-0.01 

0.26 

-0.09 

0.30 

0.00 

0.29 

0.26 

0.93 

1.72 

0.89 

106 

28,499 

0.503 

172 

0.03 

0.54 

0.26 

0.47 

-0.06 

0.52 

0.20 

0.50 

-0.55 

1.02 

1.29 

0.94 

300 

24,032 

0.430 

443 

-0.07 

0.63 

0.34 

0.71 

-0.13 

0.63 

0.26 

0.72 

0.55 

1.00 

1.31 

0.98 

OTH 

2,131,576 

43.16 

326,744 

-0.03 

0.85 

-0.03 

0.74 

-0.14 

0.86 

-0.06 

0.76 

0.19 

0.99 

0.52 

0.98 
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Table  7 

Attribution  of  unknown  699  and  644  observations. 


Catalog 

Code 

Mean  RA  (arcsec) 

Sd  RA  (arcsec) 

Mean  DEC  (arcsec) 

Sd  DEC  (arcsec) 

USNO  A2.0 

ALL 

0.0636 

0.6367 

0.2895 

0.6040 

699 

0.0243 

0.6464 

0.3470 

0.5528 

644 

0.0292 

0.3387 

0.2947 

0.3428 

UCAC2 

ALL 

-0.128 

0.5378 

0.0279 

0.5178 

USNO  B1.0 

ALL 

-0.062 

0.4869 

0.1799 

0.4434 

699 

0.0176 

0.6120 

0.2318 

0.5565 

644 

-0.052 

0.2579 

0.1839 

0.2562 

USNO  A1.0 

ALL 

-0.063 

0.7644 

-0.028 

0.7577 

Tycho-2 

ALL 

-0.006 

0.2647 

-0.001 

0.2840 

Unknown 

ALL 

-0.070 

0.7263 

0.1296 

0.6961 

699 

0.0290 

0.6704 

0.2718 

0.5995 

644 

-0.035 

0.4898 

0.3967 

0.5215 

608 

-0.085 

0.7662 

0.1077 

0.8028 

106 

-0.060 

0.5221 

0.2005 

0.5025 

689 

-0.093 

0.3019 

-0.003 

0.2915 

and  689.  The  unknown  699  residuals  and  standard  deviations  are 
an  excellent  match  for  those  of  the  A2.0  catalog;  we  feel  it  is  there¬ 
fore  extremely  likely  that  these  observations  were  reduced  using 
A2.0.  Similarly,  the  unknown  644,  608,  and  106  statistics  are  most 
consistent  with  A2.0. 

As  Table  6  indicates,  most  astrometry  reported  by  689  is  based 
on  the  Tycho-2  star  catalog;  Table  7  shows  that  the  “unknown”  689 
astrometry  is  perfectly  consistent  with  that  catalog,  making  it 
likely  these  observations  were  also  reduced  using  Tycho-2. 


5.3.  Correlation  of  residuals  with  time 

As  noted  in  Section  2,  the  statistical  error  model  proposed  by 
Carpino  et  al.  (2003)  requires  that  the  residuals  of  observations 
of  the  same  asteroid  made  by  the  same  observatory  have  only 
small  correlations  that  rapidly  decay  to  zero.  As  Fig.  11  demon¬ 
strates,  the  correlations  of  the  raw  residuals  were  not  only  signifi¬ 
cant,  but  also  persisted  for  several  days.  In  most  cases  (e.g., 
observatories  G96,  691,  and  704),  by  contrast,  the  correlations  of 
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Fig.  11.  Correlations  between  closely-spaced  observations  of  the  same  asteroid  for  several  large  observatories,  illustrating  the  reduction  in  correlation  due  to  debiasing. 
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the  debiased  residuals  were  consistently  smaller,  and  decayed  to 
insignificance  within  a  day. 

However,  the  debiased  correlations  for  observatories  703  and 
El  2  were  not  entirely  satisfactory;  although  the  debiased  correla¬ 
tions  were  indeed  smaller,  they  did  not  decay  exponentially  to 
zero.  And  as  mentioned  above,  Table  6  shows  the  debiased  RA 
and  DEC  mean  residuals  for  703  were  still  significant,  as  were 
the  debiased  DEC  mean  residuals  for  El  2.  This  led  us  to  believe  an¬ 
other  error  source  existed,  unique  to  these  observatories,  which 
might  also  account  for  the  persistent  correlations. 


5.4.  Correlation  of  residuals  with  focal  plane  position 

Observatories  703  and  El  2  both  use  the  UCAC  catalog  to  reduce 
observations;  indeed,  the  UCAC  astrometry  set  is  dominated  by 
these  two  observatories.  Referring  to  Fig.  9,  the  debiased  UCAC 
DEC  residuals  show  horizontal  striping,  while  the  debiased  UCAC 
RA  residuals  show  vertical  striping.  We  therefore  suspected  that 
there  were  significant  errors  at  the  edges  of  the  image  field,  possi¬ 
bly  caused  by  the  failure  to  account  for  high-order  field  distortions 
in  the  image  reduction  process. 

To  confirm  this,  we  used  sky  coverage  data  provided  by  obser¬ 
vatories  703  and  E12  to  deduce  where  each  asteroid  observation 
lay  in  the  image  plane.  We  then  divided  the  image  plane  into  a 
40  x  40  grid;  and  for  each  individual  pixel,  we  determined  the 


mean  debiased  RA  and  DEC  residual  for  all  observations  lying  at 
that  point  in  the  image  frame. 

As  evident  in  Fig.  12,  the  debiased  residuals  for  703  and  El  2  are 
generally  small  near  the  plate  center;  but  they  increase  rapidly 
near  the  plate  boundaries.  Similar  effects  were  also  observed  for 
observatory  G96,  which  also  uses  the  UCAC  catalog. 

Further  analysis  indicated  that  the  field  distortion  errors  for 
observatory  703  became  pronounced  between  1998  and  2003;  this 
is  consistent  with  the  introduction  of  new  optics  in  1999.  Since 
2003,  the  plate  residuals  have  been  largely  stable. 

In  Table  6,  we  note  the  significant  negative  value  of  a  for  code 
703  (Mt.  Bigelow)  with  the  UCAC  catalog,  and  the  similar  positive 
value  of  d  for  E12  (Siding  Spring).  This  may  be  related  to  the  fact 
that  the  two  installations  use  a  similar  setup,  except  that  the  cam¬ 
era  is  rotated  by  90°  at  El  2  compared  to  the  orientation  used  at 
703  (R.  McNaught,  private  communication).  Again,  the  UCAC 
astrometry  set  is  dominated  by  these  two  observatories;  and  this 
persistent  effect  is  evident  in  the  numbered  asteroid  residual  plots 
for  UCAC  (Fig.  9),  where  we  can  see  that  a  negative  residual  trend 
persists  in  the  northern  hemisphere  RA  map,  and  a  positive  signal 
is  present  in  the  southern  hemisphere  DEC  map. 

6.  Prediction  tests 

While  the  dramatic  improvement  seen  in  the  astrometric  resid¬ 
uals  when  using  the  catalog  debiasing  technique  is  encouraging,  it 
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Fig.  12.  Mean  residuals  for  the  image  planes  of  observatories  703,  El 2,  and  G96,  illustrating  the  effects  of  high  order  distortions  in  the  image  reduction  process. 
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remains  to  be  seen  whether  the  new  orbits  are  actually  superior  to 
those  obtained  with  raw  astrometry.  Adjusting  astrometry  with 
the  aim  of  reducing  the  postfit  residuals  can  be  self-defeating, 
leaving  one  with  an  orbit  that  appears  to  be  a  good  fit  but  does 
not  represent  reality.  The  deciding  test  is  whether  or  not  the  new 
orbits  actually  predict  better  or  not.  To  measure  prediction  perfor¬ 
mance,  we  look  at  how  a  series  of  orbital  fits  to  subsets  of  the  en¬ 
tire  observational  data  set  compares  to  a  “truth”  orbit  obtained  by 
fitting  over  a  very  long  observational  arc  including  all  available 
data. 

To  compare  the  prediction  performance  of  the  different  orbits 
obtained  with  raw  and  debiased  astrometry,  we  have  carefully  se¬ 
lected  a  set  of  main  belt  asteroids  with  a  thorough  CCD  observation 
history.  Specifically,  we  sought  asteroids  with  seven  consecutive 
oppositions,  each  of  which  was  well-observed  with  CCD-based 
astrometry.  To  start,  we  selected  the  5018  MBAs  with  absolute 
magnitude  H  in  the  range  14-16  and  orbital  period  in  the  range 
4-6  years.  This  absolute  magnitude  translates  to  diameters  from 
a  few  to  several  kilometers,  and  leads  to  opposition  apparent  mag¬ 
nitudes  mostly  in  the  range  17-20,  faint  enough  to  prevent  CCD 
saturation  for  most  observing  setups.  The  period  limit  enforces 
an  Earth  synodic  period  between  1.2  and  1.33  years,  ensuring  that 
there  will  have  been  about  7  oppositions  in  the  previous  10  years. 

From  among  the  initial  5018  MBAs,  we  selected  the  222  that 
actually  had  observations  at  seven  consecutive  oppositions  and 
each  opposition  with  at  least  50  CCD  observations  over  at  least  a 
50  day  time  period.  These  oppositions  generally  covered  the 
2000-2008  time  frame.  CCD  astrometry  for  which  we  could  not 
determine  the  catalog,  as  described  in  Section  3  was  not  consid¬ 
ered  in  the  opposition  astrometry  tallies. 

For  each  of  these  222  objects  we  computed  24  different  orbits, 
that  is  12  different  observational  time  spans,  each  with  raw  and 
debiased  astrometry.  All  observations  were  assigned  equal  data 
weights  of  1  arcsec.  The  first  of  the  12  fit  spans,  which  we  use  as 
proxy  for  the  truth  orbit,  was  a  fit  to  all  available  data,  including 
optical  and  CCD  observations  with  unknown  catalog.  These  cases, 
labeled  “ALL,"  used  8-12  oppositions  total,  with  median  and  mode 
of  nine  oppositions.  The  other  11  fits  were  based  on  subsets  of  the 
seven  consecutive  oppositions  described  above.  These  fits  were  in¬ 
tended  to  reveal  the  effect  of  the  debiasing  technique,  and  so  they 
used  only  CCD  observations  tied  to  a  known  catalog.  These  1 1  fits 
are  shown  schematically  in  Fig.  13,  where  each  is  shown  with  its 
identifying  tag  that  describes  the  data  set  used  in  the  respective 
fit.  As  an  example,  for  tag  “2AS"  the  “2”  indicates  that  there  are 
two  oppositions  in  the  fit,  and  the  "A”  indicates  that  the  two  oppo¬ 
sitions  are  after  the  central  (fourth)  opposition,  which  included  the 
prediction  epoch.  This  implies  that  either  oppositions  5  and  6  or  6 
and  7  were  used  in  the  fit.  The  final  character  clarifies  this  ambigu¬ 
ity:  the  “S”  in  the  tag  indicates  that  the  prediction  interval  from  the 


Table  8 

Comparison  of  relative  prediction  errors  for  222  MBAs. 


Case 

Relative  errors  (mean  ±  std.  dev.) 

Radial 

Transverse 

Normal 

1BL 

Deb. 

0.08  ±1.09 

0.22  ±1.07 

0.24  ±1.16 

Raw 

0.08  ±1.11 

0.27  ±1.09 

0.76  ±1.59 

IBS 

Deb. 

-0.22  ±1.10 

0.42  ±1.00 

0.03  ±1.07 

Raw 

-0.27  ±1.17 

0.48  ±1.05 

0.25  ±1.33 

IAS 

Deb. 

0.25  ±  0.73 

0.33  ±  0.68 

-0.1 1  ±  0.93 

Raw 

0.28  ±0.81 

0.38  ±0.77 

-0.37  ±1.64 

1AL 

Deb. 

0.06  ±  0.86 

0.38  ±  0.83 

0.23  ±  0.80 

Raw 

0.09  ±  0.94 

0.46  ±  0.91 

0.69  ±1.17 

2BL 

Deb. 

-0.01  ±1.00 

0.10  ±0.90 

0.26  ±1.09 

Raw 

0.02  ±1.13 

0.06  ±1.05 

1.03  ±2.50 

2BS 

Deb. 

-0.06  ±1.12 

0.16  ±1.18 

0.52  ±1.24 

Raw 

0.01  ±1.19 

0.09  ±1.25 

2.15  ±2.48 

2AS 

Deb. 

0.06  ±  0.85 

0.12  ±0.82 

0.36  ±1.27 

Raw 

0.13  ±0.92 

0.25  ±  0.91 

1.07  ±2.76 

2AL 

Deb. 

0.05  ±  0.83 

0.02  ±  0.80 

0.39  ±1.01 

Raw 

0.13  ±0.94 

0.09  ±  0.90 

1.20  ±1.79 

3A 

Deb. 

-0.04  ±1.12 

0.10  ±1.08 

0.47  ±1.07 

Raw 

0.07  ±1.27 

-0.01  ±1.25 

1.84  ±1.87 

3B 

Deb. 

-0.01  ±  0.82 

0.04  ±  0.78 

0.37  ±  0.93 

Raw 

0.05  ±  0.95 

0.15  ±0.91 

1.10  ±  1.75 

7 

Deb. 

-0.18  ±0.62 

-0.07  ±  0.66 

-0.02  ±  0.67 

Raw 

-0.40  ±  0.72 

-0.12  ±0.77 

0.26  ±  0.84 

fit  span  is  short,  i.e.  oppositions  5  and  6,  as  opposed  to  long  (“L”), 
which  would  denote  oppositions  6  and  7. 

We  fit  orbits  to  the  12  data  sets,  using  the  raw  astrometry  and 
then  also  debiased  astrometry.  Outlier  rejection  was  accomplished 
on  only  two  fits,  and  utilized  the  method  described  by  Carpino 
et  al.  (2003)  with  a  =  9.  First,  the  fit  for  the  debiased  ALL  case 
allowed  for  outlier  rejection  after  which  the  raw  ALL  fit  was  ob¬ 
tained  with  the  same  data  set.  Next,  the  debiased  case  7  fit  was  ob¬ 
tained  while  allowing  for  rejection,  after  which  the  outlier  set  was 
frozen  for  the  21  subsequent  fits  for  that  object.  Again,  all  data 
weights  were  1  arcsec  at  this  point. 

For  the  prediction  tests,  we  propagated  each  of  these  24  orbital 
fits,  for  each  of  the  222  asteroids,  to  the  midpoint  of  the  particular 
asteroid’s  fourth  opposition,  where  we  mapped  the  cartesian  posi¬ 
tion  and  covariance  into  the  classical  radial-transverse-normal 
(RTN)  frame  positions.  (For  consistency,  we  used  the  debiased  case 
7  RTN  rotation  for  all  the  mappings  of  a  given  asteroid.)  The  differ¬ 
ences  between  the  orbital  predictions  for  the  various  cases  and 
that  of  the  "truth”  orbit  (ALL)  in  Tables  8  and  9  allow  us  to  deter¬ 
mine  which  statistical  treatment  provided  the  better  prediction.  Of 
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Fig.  13.  Schematic  diagram  depicting  the  specific  oppositions  included  in  the  fit  spans  for  the  various  case  labels  used  in  the  prediction  tests. 
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Table  9 

Comparison  of  absolute  prediction  errors  for  222  MBAs. 


Case 

Absolute  errors  (mean  ±  std.  dev.  in  kilometers) 

Radial 

Transverse 

Normal 

1BL 

Deb. 

1274.48  ±29176.26 

19614.95  ±  115110.68 

-2.78  ±1590.12 

Raw 

728.06  ±  30500.84 

26221.85  ±120317.14 

231.44  ±1618.86 

IBS 

Deb. 

-2733.96  ±61538.02 

22526.14  ±176783.72 

-9.79  ±2152.71 

Raw 

-5056.61  ±  58828.86 

26772.96  ±162248.03 

42.21  ±  2292.45 

IAS 

Deb. 

3673.48  ±11653.44 

7868.82  ±  19486.13 

39.59  ±  687.46 

Raw 

4129.94  ±12459.56 

8323.62  ±20786.16 

-16.93  ±820.24 

1AL 

Deb. 

890.40  ±13988.90 

19178.37  ±51694.49 

72.05  ±  688.49 

Raw 

1397.80  ±15375.18 

23802.24  ±  55250.83 

198.26  ±801.11 

2BL 

Deb. 

-489.47  ±  6720.73 

3417.92  ±25080.41 

43.88  ±  256.74 

Raw 

-220.17  ±7197.23 

2199.92  ±28812.25 

189.90  ±520.15 

2BS 

Deb. 

-943.50  ±7441.94 

3754.07  ±  23252.53 

77.30  ±  280.62 

Raw 

-498.36  ±  7766.07 

2797.33  ±  23545.00 

376.65  ±527.21 

2AS 

Deb. 

313.42  ±3234.76 

757.86  ±  5524.05 

52.67  ±  220.47 

Raw 

618.90  ±3526.39 

1405.35  ±6085.32 

159.60  ±465.30 

2AL 

Deb. 

100.25  ±3887.80 

1013.54  ±  14496.75 

77.02  ±  223.32 

Raw 

438.35  ±4510.58 

2120.63  ±16640.50 

227.66  ±415.58 

3A 

Deb. 

73.93  ±1577.52 

103.78  ±3456.42 

54.49  ±134.52 

Raw 

183.94  ±1716.17 

391.43  ±3707.67 

156.76  ±  251.70 

3B 

Deb. 

-223.48  ±3921.39 

1004.72  ±10900.68 

70.24  ±183.45 

Raw 

-15.55  ±4181.81 

115.73  ±11642.35 

291.84  ±310.49 

7 

Deb. 

-9.39  ±32.78 

-8.36  ±85.96 

-3.08  ±  60.64 

Raw 

-20.25  ±37.33 

-14.26  ±  100.02 

22.16  ±  76.15 

course,  the  prediction  errors  are  statistical,  and  so  the  poorer  sta¬ 
tistical  treatment  may  provide  a  superior  prediction  for  certain  ob¬ 
jects.  The  prediction  performance  for  the  ensemble  of  222  main 
belt  asteroids  is  key  to  reaching  a  conclusion. 

As  a  typical  example,  Fig.  14  shows  the  comparison  of  the  rela¬ 
tive  prediction  errors  for  fit  2BL.  The  relative  error  is  the  error  in 
sigmas,  as  given  by  the  2BL  covariance  matrix.  In  each  plot  in 
Fig.  14  there  are  four  curves  that  compare  the  raw  and  debiased 
2BL  predictions  to  the  raw  and  debiased  ALL  orbits.  The  first  three 
plots  in  the  figure  show  the  1-D  comparison  for  each  of  the  RTN 
directions,  while  the  fourth  shows  the  3-D  position  error.  For 


reference,  the  classical  %2  cumulative  distribution  function,  which 
is  expected  for  an  uncorrelated  and  unbiased  Gaussian  error  distri¬ 
bution,  is  also  plotted. 

The  results  show  unambiguously  that  the  debiased  astrometry 
provides  predictions  that  are  both  more  accurate  and  more  in  line 
with  the  prediction  covariance.  The  greatest  improvement  is  seen 
in  the  out-of-plane  prediction  errors,  but  there  is  also  a  consistent, 
albeit  slight,  gain  in  the  radial  and  transverse  directions.  There  was 
a  possibility  of  an  ambiguous  outcome,  namely  that  each  statistical 
method  would  predict  well  against  its  own  truth  orbit  obtained,  so 
that  raw  would  compare  well  against  a  raw  truth  and  debiased 


Radial 


Transverse 


Fig.  14.  Case  2BL  cumulative  distribution  of  prediction  errors  at  the  midpoint  of  the  fourth  opposition  for  222  MBAs  described  in  the  text.  Raw  and  debiased  predictions  are 
compared  to  raw  and  debiased  truth  orbits  obtained  using  all  available  astrometry.  Abscissae  are  formal  standard  deviations,  i.e.,  sigmas. 
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would  compare  well  against  a  debiased  truth.  As  it  turns  out,  and 
as  is  clear  from  Fig.  14,  the  debiased  orbits  predict  significantly 
better  than  the  raw  orbits  no  matter  which  reference  orbit  is  used, 
and  this  is  true  for  all  1 1  cases  considered  here.  As  one  would  ex¬ 
pect,  there  is  a  slight  improvement  when  comparing  the  debiased 
orbit  to  the  debiased  truth  as  opposed  to  the  raw  truth,  but  the 
salient  point  is  that  the  debiased  orbit  outperforms  the  raw  orbit 
even  when  both  are  referenced  to  the  raw  truth  prediction. 


The  3-D  comparisons  for  all  11  fit  spans  are  shown  in  Fig.  15, 
where  it  is  clear  that  the  conclusions  from  the  2BL  case  carry  for 
every  fit  span  considered.  It  is  important  to  recognize  that  these 
plots  depict  the  relative  errors  according  to  the  computed  predic¬ 
tion  uncertainty.  The  absolute  errors  depend  heavily  on  the  arc 
length  and  prediction  interval  and  thus  vary  substantially  among 
the  different  cases.  For  example,  the  prediction  errors,  when 
measured  in  kilometers,  are  far  greater  for  1 -opposition  fits  than 
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Fig.  15.  Cumulative  distribution  of  relative  position  errors  for  predictions  at  the  midpoint  of  the  fourth  opposition.  Prediction  errors  are  obtained  from  comparison  to 
debiased  “truth”  orbits  obtained  using  all  available  astrometry.  Abscissae  in  sigmas. 
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for  2-opposition  fits.  Similarly  the  long  predictions  have  higher 
absolute  errors  than  short  ones  for  a  given  data  set.  However,  for 
a  given  object,  the  data  sample  and  weighting  were  the  same  for 
the  raw  and  debiased  fits,  and  so  the  covariance  was  the  same  also. 


6.1.  Astrometric  weighting 

While  the  discussion  so  far  has  shown  emphatically  that  orbital 
predictions  and  postfit  residuals  are  consistently  improved  with 
the  use  of  the  catalog  debiasing  technique,  some  of  the  debiased 
CDF  curves  in  Fig.  15  still  fall  below  the  x2  theoretical  prediction, 
suggesting  that  unresolved  biases  may  be  still  be  present,  or— less 
likely,  we  believe— that  the  1  arcsec  weighting  is  too  low.  Inspec¬ 
tion  of  the  plots  suggests  that  this  issue  is  mostly  associated  with 
earlier  fitspans  and  hence  earlier  data.  On  the  other  hand,  there  are 
some  cases  where  the  debiased  CDF  is  above  the  theoretical  predic¬ 
tion,  indicating  that  the  predictions  are,  in  a  statistical  sense,  too 
good,  presumably  due  to  underweighting  of  the  astrometry.  This 
pattern  appears  to  be  closely  associated  to  longer  arcs  and  more  re¬ 
cent  fit  spans. 

This  points  to  the  need  for  appropriate  weighting  among  the 
different  catalogs  and  observatories,  as  opposed  to  the  equal 
weighting  used  so  far,  which  has  the  effect  of  diminishing  the  con¬ 
tribution  of  more  accurate  observations  in  favor  of  less  accurate 
observations.  Weighting  astrometry  according  to  the  observing 
program,  star  catalog,  epoch  of  observation,  etc.  is  an  obvious  place 
to  look  for  improvement  in  orbital  predictions  (Carpino  et  al., 
2003),  and  with  this  idea  we  have  tabulated  in  Table  10  the  statis¬ 
tics  of  the  postfit  residuals  for  raw  and  debiased  astrometry  for  the 
7-opposition  fits  (i.e.,  case  7). 

Among  other  things,  this  table  affords  a  comparison  of  the  dif¬ 
ferent  catalogs  and  observing  programs,  and  is  analogous  to  Ta¬ 
ble  6,  but  for  the  smaller— and  carefully  selected— set  of  222 
MBAs.  The  first  conclusion  from  the  table,  by  now  expected,  is  that 
the  debiased  statistics  are  improved,  or  show  no  significant  change 
in  the  few  cases  where  the  raw  mean  already  happened  to  be  very 
low.  Of  particular  note  is  the  consistent  reduction  in  oi,  which  had 
not  originally  been  called  out  as  a  problem.  We  also  note  that  Ta¬ 
ble  10  is  consistent  with  the  analysis  in  Sections  5.2  and  5.4. 

But,  for  the  present  purposes,  we  propose  to  explore  the  utility 
of  applying  a  rudimentary  weighting  scheme  based  on  the  data  in 
Table  10,  as  opposed  to  the  naive  weighting  scheme  used  so  far 
that  applies  uniform  weights  to  all  astrometry.  (This  usage  of 
“naive”  follows  Carpino  et  al.  (2003)  and  should  not  be  construed 
with  its  pejorative  connotation,  but  rather  in  its  literal  sense  of 
unsophisticated  and  lacking  full  information,  as  is  certainly  the 
case  when  treating  all  observatories  equally.)  In  an  idealized  situ¬ 
ation,  with  uncorrelated  and  unbiased  astrometry,  the  a  priori  RMS 
or  errors  could  be  inferred  from  the  a  posteriori  standard  devia¬ 
tions  of  residuals.  Despite  the  fact  that  our  situation  still  falls  short 
of  the  ideal,  even  after  catalog  debiasing,  we  will  take  as  weights 
the  tabulated  values  of  (a)  and  (i5)  for  the  debiased  fits,  later  look¬ 
ing  at  the  extent  to  which  this  weighting  is  optimistic.  Specifically, 
we  use  as  weights  the  standard  deviations  for  each  catalog  and 
observatory,  or  the  “OTH”  entry  for  other  observatories,  as  listed 
on  Table  10.  The  results  from  this  rudimentary  weighting  scheme 
will  be  compared  to  that  of  the  naive  weighting  scheme. 

Comparing  the  prediction  performance  of  the  weighted  and  na¬ 
ive  fits  is  not  as  simple  as  comparing  the  raw  and  debiased  fits  be¬ 
cause  the  level  of  improvement  is  rather  more  modest.  Moreover, 
the  resulting  covariance  matrices  are  significantly  different 
between  the  naive  and  weighted  fits,  and  so  differences  in  relative 
errors  may  be  dominated  by  the  change  in  covariance.  Thus  the 
comparison  approach  demonstrated  in  Fig.  14  is  inconclusive,  be¬ 
cause  the  variations  in  the  weighted  and  naive  orbits  are  compara¬ 


ble  to  the  variations  in  the  weighted  and  naive  truth  orbits,  and  so 
there  is  not  an  adequate  truth  with  which  to  make  the  comparison. 

Because  of  this,  we  take  an  alternate  approach  to  prediction  test¬ 
ing,  namely  to  quantify  how  well  the  orbit  matches  the  fourth  oppo¬ 
sition  observations,  which  were  not  a  part  of  the  associated  orbital 
fit.  As  an  example,  case  3B  was  a  fit  to  oppositions  1-3  for  each  ob¬ 
ject.  The  RMS  of  prefit  residuals  obtained  for  opposition  4  with  the 
case  3B  initial  conditions  is  a  measure  of  how  well  the  3B  fit  pre¬ 
dicted  the  fourth  opposition.  Because  of  the  different  weighting 
strategies  we  do  not  look  at  the  RMS  of  weighted  residuals,  but 
rather  we  consider  the  raw  RMS  in  arcsec,  i.e.,  of  unnormalized  resid¬ 
uals,  to  discern  which  weighting  scheme  predicts  better. 

It  turns  out  that  there  is  a  clear  but  modest  improvement  in  the 
absolute  predictions  obtained  with  this  rudimentary  weighting. 
Fig.  16  gives  the  cumulative  distributions  of  the  fourth  opposition 
RMS  of  unweighted  residuals.  Curves  are  presented  for  the  RA  and 
DEC  RMS,  combined  and  separate,  for  each  of  the  11  fitspans.  Over¬ 
all,  the  DEC  predictions  are  superior  to  the  RA  predictions,  with  the 
combined  RMS  falling  somewhere  in  between.  This  was  antici¬ 
pated  because  most  orbit  prediction  uncertainty  accumulates 
along  track,  which  maps  mostly  into  RA.  The  case  7  plot  is  included 
although  in  this  case  the  fourth  opposition  RMS  is  postfit  and  so  is 
not  a  test  of  prediction  performance.  Among  the  10  other  cases,  we 
see  that  the  weighted  distributions  are  superior  to  the  naive  distri¬ 
butions,  or  they  are  very  close,  as  is  the  case  for  case  2BS.  Thus  we 
conclude  that  the  weighting  does  produce  better  orbital  predic¬ 
tions.  This  conclusion  is  reinforced  quantitatively  by  Table  11, 
which  lists  the  median  and  90th  percentile  combined  RMS  for  each 
of  these  ten  cases.  In  every  case,  the  weighted  RMS  is  either  signif¬ 
icantly  less— or  substantially  the  same  as— the  naive  RMS.  The  table 
also  lists  the  geometric  mean  and  the  median  of  the  ratios  between 
the  weighted  and  naive  RMS  values,  where  a  number  less  than  one 
indicates  a  superior  weighted  prediction.  It  is  clear  again  from  the 
table  that  the  weighting  produces  consistently  better  predictions. 
We  note  that  the  geometric  mean  is  less  than  the  median  in  9  of 
the  10  cases,  implying  that  for  individual  objects  with  a  large  dis¬ 
parity  between  the  naive  and  weighted  predictions,  the  results 
tend  to  favor  the  weighted  solutions,  a  conclusion  supported  by 
inspection  of  the  raw  data. 

While  we  have  established  that  the  rudimentary  weighting 
scheme  does  produce  predictions  that  are  superior  to  the  naive 
scheme,  at  least  in  an  absolute  sense,  what  is  so  far  missing  from 
the  analysis  is  an  understanding  of  how  well  the  prediction  errors 
match  the  a  posteriori  covariance.  As  one  would  expect,  unresolved 
biases  and  correlations  cause  the  standard  deviations  from  Table  1 1 
to  be  significantly  optimistic  when  translated  into  weights,  and  the 
resulting  a  posteriori  orbit  covariances  significantly  underestimate 
the  orbit  estimation  errors.  To  accommodate  this  we  compute  and 
tabulate  the  safety  factor  needed  to  scale  the  covariance  and  bring 
it  into  line  with  the  orbit  errors  for  each  case. 

As  detailed  in  Fig.  17,  with  case  3A  used  as  an  example,  this 
safety  factor  was  obtained  from  the  ratio  of  the  normalized  posi¬ 
tion  errors  to  that  expected  from  an  idealized  Gaussian  distribu¬ 
tion.  The  values  listed  in  Table  1 1  were  taken  at  the  median,  but 
as  is  clear  from  the  figure,  the  value  is  reasonably  constant  across 
the  distribution.  Six  of  the  values  fall  within  or  close  to  the  range 
1.5-1. 8,  while  the  other  four  are  larger,  roughly  in  the  range  2.3- 
3.0.  It  may  be  significant  that  these  four  cases  with  large  safety  fac¬ 
tors  match  exactly  the  four  cases  that  include  data  from  opposition 
2  (see  Fig.  13).  In  rough  terms,  this  corresponds  to  the  period  from 
late  2001  to  late  2002,  which  happens  to  be  a  time  frame  when  the 
search  efforts  were  particularly  dominated  by  code  704,  the  LIN¬ 
EAR  search  program  (Stokes  et  al.,  2000),  with  approximately 
80%  of  astrometry  obtained  by  LINEAR.  It  is  also  the  time  frame 
when  the  program’s  second  telescope  was  being  brought  into  oper¬ 
ation,  which,  incidentally,  we  are  not  able  to  distinguish  from  their 


Table  10 

Residual  statistics  for  222  main  belt  asteroids  by  catalog  and  observatory. 


Code 

Use/delete 

Percent 

Debiased  (arcsec) 

Raw  (arcsec) 

Ratio  debiased/raw 

a 

(a) 

3 

<<5> 

a 

(a> 

3 

(S) 

a 

(a) 

3 

w 

USNO  A2.0 

ALL 

103,471/440 

68.4 

0.03 

0.60 

0.10 

0.58 

0.06 

0.62 

0.34 

0.60 

0.45 

0.96 

0.29 

0.97 

704 

85,007/380 

82.2 

0.04 

0.62 

0.11 

0.60 

0.08 

0.64 

0.34 

0.62 

0.56 

0.96 

0.31 

0.97 

699 

4343/6 

4.2 

-0.01 

0.46 

0.06 

0.39 

0.03 

0.49 

0.31 

0.41 

-0.20 

0.95 

0.19 

0.96 

691 

3717/11 

3.6 

-0.06 

0.31 

0.01 

0.34 

0.00 

0.33 

0.25 

0.38 

11.99 

0.95 

0.05 

0.89 

608 

3542/13 

3.4 

0.06 

0.63 

0.09 

0.76 

0.08 

0.67 

0.39 

0.78 

0.69 

0.94 

0.24 

0.97 

703 

3343/22 

3.2 

-0.17 

0.62 

0.01 

0.57 

-0.12 

0.64 

0.30 

0.58 

1.50 

0.96 

0.02 

0.97 

644 

1641/0 

1.6 

-0.03 

0.24 

0.11 

0.28 

0.00 

0.31 

0.37 

0.30 

19.67 

0.77 

0.30 

0.92 

OTH 

1878/8 

1.8 

-0.02 

0.51 

-0.01 

0.40 

0.02 

0.52 

0.24 

0.41 

-1.12 

0.97 

-0.04 

0.95 

UCAC1  &  2 

ALL 

27,211/70 

18.0 

-0.12 

0.43 

0.07 

0.42 

-0.17 

0.44 

0.05 

0.45 

0.70 

0.98 

1.27 

0.94 

703 

12,978/37 

47.7 

-0.21 

0.49 

0.09 

0.46 

-0.26 

0.50 

0.06 

0.47 

0.80 

0.98 

1.62 

0.96 

G96 

6787/5 

24.9 

-0.04 

0.25 

-0.05 

0.21 

-0.09 

0.27 

-0.10 

0.25 

0.46 

0.92 

0.47 

0.84 

E12 

4891/13 

18.0 

-0.03 

0.41 

0.17 

0.43 

-0.07 

0.42 

0.27 

0.44 

0.36 

0.99 

0.62 

0.98 

683 

1033/13 

3.8 

0.00 

0.60 

0.09 

0.78 

-0.06 

0.62 

0.06 

0.78 

-0.02 

0.98 

1.41 

0.99 

OTH 

1522/2 

5.6 

-0.02 

0.33 

0.01 

0.30 

-0.06 

0.34 

-0.02 

0.33 

0.39 

0.97 

-0.21 

0.89 

USNO  B1.0 

ALL 

16,056/19 

10.6 

0.02 

0.37 

0.10 

0.35 

-0.08 

0.37 

0.19 

0.37 

-0.27 

0.98 

0.50 

0.93 

699 

9341/8 

58.2 

0.09 

0.42 

0.14 

0.40 

-0.01 

0.42 

0.25 

0.42 

-8.31 

0.99 

0.56 

0.96 

644 

4092/0 

25.5 

-0.02 

0.18 

0.04 

0.17 

-0.13 

0.20 

0.14 

0.23 

0.18 

0.90 

0.30 

0.71 

691 

2165/6 

13.5 

-0.20 

0.25 

0.02 

0.28 

-0.30 

0.26 

0.07 

0.31 

0.66 

0.94 

0.27 

0.91 

OTH 

458/5 

2.9 

-0.01 

0.49 

0.04 

0.40 

-0.10 

0.50 

0.15 

0.43 

0.05 

0.99 

0.29 

0.94 

GSC-ACT 

ALL 

1628/19 

1.1 

0.06 

0.55 

0.12 

0.56 

-0.04 

0.55 

0.10 

0.58 

-1.52 

1.01 

1.15 

0.97 

333 

1023/5 

62.8 

0.09 

0.52 

0.06 

0.46 

-0.05 

0.52 

0.06 

0.49 

-1.87 

1.00 

1.07 

0.93 

OTH 

605/14 

37.2 

0.02 

0.60 

0.21 

0.70 

-0.03 

0.59 

0.18 

0.71 

-0.53 

1.01 

1.20 

0.99 

USNOAl.O 

ALL 

1315/10 

0.9 

-0.02 

0.58 

0.01 

0.51 

-0.11 

0.62 

-0.05 

0.57 

0.21 

0.94 

-0.20 

0.89 

704 

1103/10 

83.9 

0.00 

0.59 

0.00 

0.52 

-0.08 

0.63 

-0.06 

0.59 

0.00 

0.93 

0.07 

0.88 

OTH 

212/0 

16.1 

-0.14 

0.55 

0.09 

0.47 

-0.25 

0.55 

-0.02 

0.46 

0.59 

0.99 

-4.26 

1.03 

CMC 

ALL 

974/0 

0.6 

0.02 

0.44 

0.04 

0.36 

-0.03 

0.44 

-0.02 

0.37 

-0.73 

0.99 

-1.73 

0.96 

Tycho-2 

ALL 

459/5 

0.3 

0.04 

0.40 

-0.04 

0.40 

-0.04 

0.42 

-0.01 

0.42 

-0.81 

0.95 

2.47 

0.94 

689 

370/0 

80.6 

0.01 

0.26 

-0.04 

0.31 

-0.07 

0.28 

-0.01 

0.36 

-0.20 

0.91 

4.38 

0.88 

OTH 

89/5 

19.4 

0.12 

0.72 

-0.04 

0.63 

0.05 

0.75 

-0.04 

0.62 

2.32 

0.97 

0.98 

1.02 

GSC  1 

ALL 

132/1 

0.1 

-0.01 

0.45 

0.11 

0.44 

-0.04 

0.45 

0.05 

0.45 

0.32 

0.99 

2.18 

0.97 

Unknown 

ALL 

9515 

N/A 

0.07 

0.58 

0.34 

0.59 

-0.03 

0.57 

0.31 

0.59 

-2.46 

1.03 

1.13 

0.99 

699 

3090 

32.5 

0.05 

0.56 

0.44 

0.41 

-0.04 

0.53 

0.40 

0.41 

-1.23 

1.04 

1.10 

0.99 

644 

2645 

27.8 

0.09 

0.43 

0.41 

0.49 

-0.01 

0.40 

0.36 

0.50 

-8.47 

1.08 

1.13 

0.98 

608 

2049 

21.5 

0.13 

0.77 

0.31 

0.85 

0.02 

0.76 

0.29 

0.84 

7.28 

1.01 

1.05 

1.00 

D29 

484 

5.1 

0.03 

0.36 

-0.05 

0.32 

-0.03 

0.37 

-0.10 

0.34 

-1.13 

0.97 

0.46 

0.95 

689 

206 

2.2 

-0.02 

0.23 

-0.01 

0.29 

-0.07 

0.26 

-0.05 

0.33 

0.32 

0.90 

0.18 

0.86 

106 

178 

1.9 

0.03 

0.56 

0.33 

0.51 

-0.06 

0.55 

0.28 

0.52 

-0.55 

1.02 

1.18 

0.98 

300 

177 

1.9 

-0.07 

0.60 

0.28 

0.69 

-0.14 

0.61 

0.24 

0.71 

0.51 

1.00 

1.18 

0.97 
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Case  1 BL 


Case  2BL 


Case  3B 
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Case  2BS 


Case  3A 


Case  IAS 


Case  1  AL 


Case  2AS 
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Case  7 


- Naive 

- Weighted 

-  -  -  Naive,  RA 

- Weighted,  RA 

- Naive,  DEC 

—  Weighted,  DEC 


Fig.  16.  Cumulative  distribution  of  prefit  RMS  of  fourth  opposition  residuals  for  various  orbital  solutions  with  222  MBAs.  The  comparison  is  between  the  naive  data  weighting 
scheme  and  a  rudimentary  scheme  based  on  observatory  and  catalog.  For  each  plot,  the  absicca  is  the  prefit  RMS  in  arcsec  and  the  ordinate  is  the  cumulative  distribution.  Case 
7  is  included  for  reference,  although  only  postfit  residuals  are  available  for  the  fourth  opposition. 


first  telescope  by  reference  to  the  astrometric  data  delivered  by  the 
MPC.  In  any  case,  the  mean  of  the  tabulated  safety  factors  is 
2.03  ±  0.56.  From  all  of  this  we  conclude  that  scaling  the  a  posteri¬ 
ori  standard  deviation  of  residuals  by  a  factor  of  2,  and  using  these 
as  weights  leads  to  improved  absolute  and  relative  predictions,  at 
least  for  the  present  ensemble  of  222  MBAs. 


7.  Applications 

7.1.  Asteroid  mass  estimation 

Careful  asteroid  orbit  determination  affords  the  possibility  of 
observing  the  perturbation  by  many  of  the  larger  main  belt 
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Table  11 

Comparison  of  prefit  RMS  of  fourth  opposition  residuals,  with  and  without  observatory-specific  weighting,  for  222  MBAs. 


Case 

Ratio  weighted/naive 

Median  (arcsec) 

90th  percentile  (arcsec) 

Safety  factor 

Geom.  mean 

Median 

Naive 

Weighted 

Naive 

Weighted 

1AL 

0.715 

0.729 

14.979 

10.014 

48.179 

34.651 

1.51 

IAS 

0.835 

0.839 

6.181 

5.076 

18.042 

15.558 

1.49 

IBS 

0.884 

0.949 

13.591 

12.597 

49.820 

52.295 

1.82 

1BL 

0.834 

0.890 

27.926 

23.061 

80.177 

74.104 

2.29 

2AL 

0.726 

0.736 

4.748 

3.403 

12.275 

8.495 

1.59 

2AS 

0.780 

0.818 

1.928 

1.379 

5.367 

3.589 

1.79 

2BS 

0.936 

0.947 

3.557 

3.565 

10.168 

9.112 

2.97 

2BL 

1.008 

0.964 

6.186 

5.600 

18.514 

16.634 

2.65 

3A 

0.862 

0.920 

1.207 

1.059 

3.227 

2.644 

1.54 

3B 

0.936 

0.969 

2.324 

2.026 

5.637 

5.827 

2.66 

Normalized  Position  Error  and  Ratio 

Fig.  17.  Cumulative  distribution  of  normalized  position  errors  for  weighted  fits  to 
case  3A  (as  an  example),  depicted  with  the  A/2  CDF.  The  ratio  between  the  actual 
and  ideal  errors  at  a  given  cumulative  distribution  level  is  also  depicted.  This  is  the 
scale  factor  required  to  bring  the  weighted  a  posteriori  covariance  into  agreement 
with  the  Gaussian  distribution. 


asteroids  on  some  of  the  more  numerous  smaller  asteroids,  termed 
test  particles  in  this  context.  The  combination  of  the  measured  per¬ 
turbation,  usually  an  offset  in  orbital  longitude,  and  the  estimated 
approach  circumstances  leads  to  an  estimate  of  the  mass  of  the  lar¬ 
ger  asteroid.  However,  the  mass  estimate  can  be  skewed  if  the 
asteroid  orbit  determinations  are  subject  to  catalog  bias. 

Aside  from  events  associated  with  1  Ceres  or  4  Vesta,  most 
asteroid  mass  determination  events  thus  far  discovered  involve 
cumulative  deflections  on  the  order  of  a  few  arcsec  (or  less).  Since 
the  debiasing  corrections  are  non-trivial  on  this  scale,  we  expected 
some  improvement  in  mass  determination  relative  to  the  use  of 
raw  observations. 

Having  said  that,  successful  mass  determination  encounters 
usually  involve  changes  in  an  asteroid’s  semimajor  axis  (and  thus 
its  mean  motion),  resulting  in  cumulative  along-track  shifts  in 
sky  position  versus  the  pre-encounter  trajectory.  As  we  saw  in 
Fig.  7,  however,  while  debiasing  does  indeed  improve  determina¬ 
tion  of  the  semimajor  axis,  the  improvement  is  somewhat  less  pro¬ 
nounced  than  for  the  eccentricity  or  inclination.  Therefore,  it 
seemed  likely  that  the  improvement  in  mass  determination  would 
be  relatively  modest. 

Table  12  illustrates  the  results  of  1 1  such  comparative  trials,  rep¬ 
resenting  all  known  mass  determination  events  (aside  from  1  Ceres 
or  4  Vesta)  with  encounters  after  1997,  and  with  a  significance  (de¬ 
fined  as  the  ratio  of  mass  to  mass  uncertainty)  greater  than  2.5.  The 
“best”  mass  values  (representing  published  values)  were  calculated 


using  the  full  set  of  available  observations;  while  the  “debiased” 
and  “raw”  masses  were  calculated  using  only  post-1995  observa¬ 
tions  (since  observations  before  this  time  precede  availability  of 
the  USNO  catalogs,  and  thus  cannot  be  debiased).  Consequently, 
the  “best”  masses  benefit  from  a  far  better  knowledge  of  the  test 
asteroid’s  pre-encounter  trajectory.  However,  our  intent  here  is 
not  to  obtain  the  best  possible  mass  estimates,  but  rather  to  deter¬ 
mine  if  debiasing  results  in  superior  values. 

We  frankly  believe  that  the  formal  uncertainties  in  asteroid 
mass  determination  likely  understate  the  actual  effects  of  unmod¬ 
eled  perturbations  and  unknown  error  sources;  so  while  we  have 
listed  the  best  currently  available  mass  for  each  asteroid  (with 
no  debiasing  applied),  these  values  may  not  be  sufficiently  well 
determined  to  act  as  absolute  references  for  comparison.  We  have 
therefore  also  provided  estimates  of  bulk  density  and  porosity; 
clearly  unrealistic  values  should  act  as  an  additional  check  in  these 
comparisons. 

In  10  of  the  11  cases,  the  debiased  mass  is  closer  to  the  best  cur¬ 
rently  available  value  than  the  mass  derived  from  raw  observa¬ 
tions.  The  exception  being  the  case  of  10/24,433,  where  the  raw 
and  debiased  results  are  virtually  identical.  This  advantage  is  most 
dramatically  apparent  in  the  case  of  9/29,818,  where  the  difference 
between  the  debiased  and  raw  masses  is  on  the  order  of  ten  times 
the  formal  uncertainty  (i.e.,  10-cr).  Indeed,  one  might  argue  that  the 
debiased  mass  for  9  Metis,  indicating  a  near-zero  bulk  porosity,  is 
more  realistic  than  the  best  currently  available  value,  which  results 
in  a  negative  porosity,  which  we  note  again  is  derived  with  raw 
astrometry. 

The  apparent  advantage  of  using  debiased  observations  is  still 
evident,  although  not  as  starkly,  in  the  results  for  1/13,801,  7/ 
52,443,  10/3946,  15/50,278,  and  189/6224.  In  each  case,  not  only 
is  the  debiased  mass  closer  to  the  best  available  value,  but  the 
resulting  bulk  porosity  is  equally  (if  not  more)  reasonable.  In  the 
cases  of  3/17,715,  6/4497  and  10/57,493,  the  formal  uncertainty 
in  the  calculated  mass  is  larger  than  the  mass  itself;  nevertheless, 
we  believe  it  is  significant  that  the  calculated  masses  are  consis¬ 
tently  closer  to  the  best  available  values,  and  the  resulting  bulk 
porosities  are  more  reasonable  than  those  resulting  from  raw 
observations. 


7.2.  Impact  hazard  assessment 

From  a  public  perspective,  one  of  the  most  notable  applications 
of  asteroid  orbit  determination  is  the  assessment  of  the  impact  risk 
posed  by  potentially  hazardous  asteroids.  A  natural  question  is  to 
what  extent  the  debiasing  techniques  proposed  in  this  paper  will 
affect  such  analyses.  There  are  two  independent,  automated 
impact-monitoring  systems,  Sentry5  at  NASA’s  Jet  Propulsion 


5  http://neo.jpl.nasa.gov. 
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Table  12 

Mass/porosity  determinations  using  raw  and  debiased  observations. 


Method 

Subject  Ast. 

Radius  (km) 

Test  ast. 

Arc 

Nobs 

Enc.  yr. 

Mass  (M0) 

Rel.  Diff.a 

Bulk  density 
(g/cm)e 

Bulk  porosityb 

Raw 

1  Ceres 

476.2  ±  3.6d 

13801 

1997.8-2008.9 

451 

1998.1 

4.42  ±1.16 

X 

10-’° 

1.94  ±0.51 

0.43  ±0.15 

Debiased 

1997.8-2008.9 

451 

1998.1 

4.76  ±1.16 

X 

to-10 

0.29 

2.09  ±0.51 

0.38  ±0.15 

Bestc 

4.753  ±0.007  x  10~10 

2.09  ±  0.05 

0.39  ±  0.01 

Raw 

3  Juno 

129  ±  3e,f 

17715 

1997.9-2008.9 

497 

2000.8 

2.71  ±1.99 

X 

10-" 

6.02  ±  4.44 

-0.77  ±1.31 

Debiased 

1997.9-2008.9 

497 

2000.8 

1.46  ±1.98 

X 

to-11 

-0.63 

3.25  ±  4.40 

0.05  ±1.30 

Bestc 

1.51  ±0.03 

X 

ltr” 

3.36  ±  0.24 

0.01  ±  0.07 

Raw 

6  Hebe 

93±1.45hJ 

4497 

1997.1-2008.5 

654 

1997.5 

22.30  ±25.30  : 

<  to-12 

13.17  ±14.95 

-2.87  ±  4.40 

Debiased 

1997.1-2008.5 

654 

1997.5 

7.20  ±  25.20  : 

<  to-12 

-0.60 

4.25  ±  14.88 

-0.25  ±4.38 

Best8 

6.46  ±  0.32 

X 

to-12 

3.81  ±0.26 

-0.12  ±0.08 

Raw 

7  Iris 

106  ±  5th 

52443 

1998.8-2008.6 

396 

2001.2 

12.1  ±5.96 

X 

to-12 

4.80  ±  2.46 

-0.43  ±  0.73 

Debiased 

1998.8-2008.6 

396 

2001.2 

8.42  ±  5.90 

X 

to-12 

-0.63 

3.33  ±  2.38 

0.00  ±  0.71 

Best8 

6.86  ±  0.50 

X 

to-12 

2.72  ±  0.43 

0.19  ±0.13 

Raw 

9  Metis 

87  ±  6hJ 

29818 

1996.5-2008.9 

468 

2005.2 

13.1  ±6.33 

X 

to-12 

9.47  ±  4.98 

-1.83  ±1.49 

Debiased 

1996.5-2008.9 

468 

2005.2 

4.54  ±  0.86 

X 

to-12 

-10.0 

3.28  ±  0.92 

0.02  ±  0.27 

Best8 

5.70  ±1.40 

X 

to-12 

4.12  ±1.33 

-0.23  ±  0.40 

Raw 

10  Hygiea 

215.3  ±  3.4th 

3946 

1996.8-2008.9 

1093 

1998.2 

7.15  ±5.80 

X 

to-11 

3.40  ±  2.77 

-0.15  ±0.94 

Debiased 

1996.8-2008.9 

1093 

1998.2 

6.51  ±  5.80 

X 

10-” 

-0.11 

3.10  ±2.77 

-0.05  ±  0.94 

Best8 

4.45  ±  0.07 

X 

10-” 

2.12  ±0.11 

0.28  ±  0.04 

Raw 

10  Hygiea 

215.3  ±  3.4th 

24433 

1995.1-2008.9 

504 

2004.5 

4.11  ±0.85 

X 

10-” 

1.96  ±0.42 

0.34  ±0.14 

Debiased 

1995.1-2008.9 

504 

2004.5 

4.08  ±  0.83 

X 

10-” 

-0.04 

1.94  ±0.41 

0.34  ±0.14 

Best8 

4.45  ±  0.07 

X 

10-” 

2.12  ±0.11 

0.28  ±  0.04 

Raw 

10  Hygiea 

215.3  ±  3.4th 

57493 

1996.4-2007.5 

314 

2000.6 

2.41  ±  5.90 

X 

10-” 

1.15  ±  2.81 

0.61  ±  0.95 

Debiased 

1996.4-2007.5 

314 

2000.6 

3.80  ±7.13 

X 

10-” 

0.19 

1.81  ±3.40 

0.39  ±1.15 

Best8 

4.45  ±  0.07 

X 

10-” 

2.12  ±0.11 

0.28  ±  0.04 

Raw 

10  Hygiea 

215.3  ±3.4r,-h 

75794 

2000.0-2008.2 

208 

2005.7 

5.87  ±4.28 

X 

10-” 

2.80  ±2.04 

0.05  ±  0.69 

Debiased 

2000.0-2008.2 

208 

2005.7 

4.06  ±  3.96 

X 

10-” 

-0.46 

1.94  ±1.89 

0.34  ±  0.64 

Best8 

4.45  ±  0.07 

X 

10-” 

2.12  ±0.11 

0.28  ±  0.04 

Raw 

15  Eunomia 

134.1  ±  7.5th 

50278 

1998.9-2008.2 

220 

2002.2 

1.65  ±0.13 

X 

10-” 

3.24  ±  0.60 

0.03  ±0.18 

Debiased 

1998.9-2008.2 

220 

2002.2 

1.59  ±0.13 

X 

10-” 

-0.41 

3.14  ±0.59 

0.06  ±0.17 

Best8 

1.57  ±0.02 

X 

10-” 

3.09  ±  0.52 

0.08  ±0.16 

Raw 

189  Phthia 

19  ±  2h 

6224 

1995.0-2008.9 

547 

1999.1 

7.93  ±  2.04 

X 

10-’4 

5.49  ±  2.24 

-0.64  ±  0.67 

Debiased 

1995.0-2008.9 

547 

1999.1 

7.06  ±  2.20 

X 

10-14 

-0.40 

4.89  ±2.17 

-0.46  ±  0.65 

Best8 

1.93  ±0.41 

X 

10-’4 

1.34  ±0.51 

0.60  ±0.15 

a  The  relative  difference  is  the  difference  between  the  raw  and  debiased  mass  estimates,  normalized  by  the  stated  debiased  (formal)  uncertainty. 
b  We  assume  grain  densities  of  2.4  g/cme  for  S  types  (i.e.,  asteroids  3,  6,  7,  9,  15  and  189)  and  G  types  (i.e.,  1  Ceres),  and  3.0  g/cme  for  C  types  (i.e.,  10  Hygiea). 
c  References:Pitjeva  (2005). 
d  Thomas  et  al.  (2005). 
e  Millis  et  al.  (1981). 
f  Kaasalainen  et  al.  (2002). 
g  Baer  et  al.  (2008). 
h  Tedesco  et  al.  (2002). 

1  Torppa  et  al.  (2003). 

J  Thomas  et  al.  (1996). 


laboratory  and  NeoDys6  at  the  University  of  Pisa  (Milani  et  al.,  2005), 
which  continually  monitor  the  ever-changing  asteroid  orbit  catalog 
for  impact  possibilities  and  assess  the  associated  probabilities.  But 
to  date  these  assessments  have  been  obtained  with  raw,  and  gener¬ 
ally  unweighted,  astrometry,  and  so  these  assessments  will  need  to 
be  re-run  with  catalog  bias  corrections  applied.  While  we  do  not  ex¬ 
pect  the  new  assessments  to  wholly  invalidate  the  majority  of  cur¬ 
rent  potential  impacts,  it  will  undoubtedly  lead  to  many  impacts 
possibilities  being  eliminated  as  no  longer  consistent  with  the  new 
(debiased)  astrometry,  and  to  many  new  potential  impacts  that  were 
previously  not  compatible  with  the  original  (raw)  astrometry.  For 
many  cases  that  persist,  there  will  be  significant  corrections  to  the 
tabulated  probability  of  impact. 

7.2.3.  99942  Apophis 

Of  particular  interest  is  the  Asteroid  99942  Apophis,  discovered 
by  R.  Tucker,  D.  Tholen  and  F.  Bernardi  at  Kitt  Peak  in  June  2004, 
and  subsequently  lost  due  to  weather  and  low  solar  elongation. 
It  was  serendipitously  recovered  by  G.  Garradd  at  Siding  Spring 
in  December  2004  (Gilmore  et  al„  2004).  in  the  days  following 
the  rediscovery,  it  became  clear  that  there  would  be  a  very  close 
Earth  approach  in  April  2029,  and  at  least  initially,  the  possibility 


6  http://newton.dm.unipi.it/neodys/. 


of  impact  could  not  be  ruled  out.  The  impact  probability  reached 
a  peak  of  2.7%  (Chesley,  2006)  at  which  point  pre-discovery  mea¬ 
surements  from  March  2004  were  reported  (Gleason  et  al.,  2004), 
which  allowed  the  2029  impact  to  be  ruled  out.  Still,  the  2029  close 
approach  would  be  very  close  indeed,  nominally  at  about 
38,000  km  from  the  geocenter  on  April  13,  2029. 

At  this  point  the  key  question  became  what  would  happen  after 
2029,  when  the  orbit  would  be  deflected  significantly,  from  the 
Aten  to  the  Apollo  class,  and  the  orbital  uncertainty  would  expand 
dramatically  due  to  the  scattering  effect  of  the  Earth  close  ap¬ 
proach.  Resonant  return  impacts  (Milani  et  al„  1999)  were  possible 
in  the  years  following  2029,  but  many  of  these  were  ruled  out  as 
the  2029  b-plane  uncertainty  region  was  reduced  through  contin¬ 
uing  observations.  As  of  this  writing,  only  the  2036  and  2037  po¬ 
tential  impacts  remain  (7:6  and  8:7  resonances,  respectively). 
These  impacts  will  take  place  if  the  true  asteroid  b-plane  intersec¬ 
tion  falls  on  a  narrow  region,  called  a  keyhole,  which  can  be  de¬ 
fined  as  the  intersection  of  the  uncertainty  region  and  the 
Valsecchi  circle  associated  with  the  given  resonance  (Valsecchi 
et  al„  2003).  The  key  question  here  is  how  much  the  current  fa- 
plane  footprint  is  affected  by  the  introduction  of  catalog  debiasing 
to  the  orbit  determination. 

The  current  Apophis  orbit  determination,  JPL  solution  140,  is 
obtained  with  raw  astrometry  and  uses  manually  adjusted  data 
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Fig.  18.  For  99942  Apophis  orbital  fits,  effect  of  debiasing  on  postfit  residuals  for 
Mauna  Kea  (568)  and  Kitt  Peak  (695)  astrometry.  One  sigma  covariance  ellipses  for 
the  respective  data  sets  are  shown. 

weighting  reflective  of  rough  estimates  of  the  observation  quality. 
The  data  set  at  present  uses  738  optical  measurements  from  2004- 
03-15  to  2006-08-16,  plus  two  radar  delay  and  Doppler  measure¬ 
ments  (Giorgini  et  al.,  2008).  We  have  also  performed  a  fit  with 
debiased  astrometry,  using  the  same  data  set  and  weights  as  solu¬ 
tion  140.  The  postfit  residuals  for  the  raw  solution  show  a  signifi¬ 
cant  positive  offset  in  the  declination  residuals,  especially  for  the 
later  observations  derived  by  Tholen  et  al.  (2008).  This  trend 
apparently  continues,  as  Tholen  et  al.  (2008)  report  that  they  have 
reduced  close  to  200  additional,  and  so  far  unreported,  observa¬ 
tions  through  January  2008,  which  show  a  persistent  positive  dec¬ 
lination  bias.  It  is  clear  from  inspection  of  the  data  that  the 
astrometric  reductions  by  Tholen  et  al.  (2008),  consisting  of  the 
observations  made  at  Mauna  Kea  Observatory  (MKO,  code  568) 
as  well  as  the  discovery  series  of  measurements  from  Kitt  Peak  Na¬ 
tional  Observatory  (KPNO,  code  695),  have  an  internal  consistency 
far  superior  to  other  contributors  of  optical  astrometry. 

Through  the  MPC,  we  currently  have  access  to  a  combined  total 
of  only  33  KPNO  and  MKO  observations.  The  residuals  for  these 
observations,  for  both  the  raw  and  debiased  fits,  are  plotted  in 
Fig.  18.  The  plot  shows  the  clear  bias  in  the  astrometry  from  the 
raw  fit,  where  there  is  not  even  a  single  negative  declination  resid¬ 
ual,  and  the  mean  is  located  at  -0.5  and  +1.5  standard  deviations 
in  RA  and  DEC,  respectively.  In  contrast,  the  residuals  from  the 
debiased  fit  show  effectively  zero  mean  and  reduced  scatter  as 
indicated  by  the  smaller  covariance  ellipse.  It  is  important  to  rec¬ 
ognize  that  we  have  only  plotted  residuals  for  33  of  the  738  obser¬ 
vations  used  in  the  fit,  but  the  debiased  fit  has  applied  catalog  bias 
corrections,  according  to  the  technique  described  earlier,  to  all 
observations.  Thus  the  effect  of  debiasing  seen  in  Fig.  18  comes 
from  two  directions:  the  orbit  is  improved  so  that  orbital  error 
no  longer  contributes  to  the  declination  bias  problem,  and  the 
MKO  and  KPNO  observations  themselves  are  improved  through 
debiasing.  Comparable  plots  for  other  data  sets  are  not  as  illumi¬ 
nating  because  none  have  such  small  internal  noise  levels,  and  so 
the  associated  ellipses  are  much  larger  than  those  seen  in  Fig.  18. 
Even  so,  the  clear  conclusion,  which  is  to  be  expected  by  now,  is 
that  the  debiased  astrometry  allows  for  more  consistent  fits,  which 
should  lead  to  better  predictions. 


Fig.  19  reveals  the  effect  of  debiasing  on  the  2029  b-plane  inter¬ 
section,  and  specifically  the  keyhole  overlap.  We  use  the  b-plane 
formalism  and  notation  described  by  Valsecchi  et  al.  (2003),  where 
the  Earth’s  inertial  velocity  projects  onto  the  -f-axis.  The  raw  and 
debiased  Apophis  prediction  ellipses,  with  one-sigma  semiaxes  of 
15  x  560  km,  are  noticeably  offset  from  each  other.  The  offsets 
are  (Ai,  AC)  ~  (23,-170)  km,  which  corresponds  to  1.7er  in  {  and 
0.3ff  in  L  While  the  {  offset  is  small  in  absolute  terms,  it  is  by  far 
the  larger  in  a  relative  sense,  and  this  is  because  the  catalog  bias 
is  dominated  by  declination  errors,  and  so  leads  to  errors  in  incli¬ 
nation,  which  map  primarily  into  the  c  direction.  From  the  figure, 
we  can  see  that  the  effect  of  debiasing  on  the  2029  b-plane  predic¬ 
tion  is  to  move  the  prediction  region  away  from  the  2036  Valsecchi 
circle,  downgrading  the  threat  from  the  associated  potential  im¬ 
pact,  but  slightly  upgrading  that  from  the  2037  impact.  We  refer 
to  the  long  axis  of  the  ellipse  as  the  Line  of  Variations  (LOV),  and 
it  is  informative  to  place  the  keyholes  in  the  frame  of  the  LOV,  mea¬ 
sured  in  sigmas.  The  LOV  positions,  and  associated  impact  proba¬ 
bilities,  computed  as  described  in  Milani  et  al.  (2002)  are 
tabulated  in  Table  13. 

While  we  do  not  at  this  point  have  access  to  the  unpublished 
astrometry  described  by  Tholen  et  al.  (2008),  we  can  coarsely  esti¬ 
mate  the  effect  of  these  data  by  simulating  16  nights  of  data,  with 
three  positions  each  night,  totaling  48  positions.  The  16  nights  se¬ 
lected  are  the  15th,  16th,  17th  and  18th  of  December  2008,  January 
2007,  December  2007  and  January  2008,  and  we  weight  the  simu¬ 
lated  data  at  0.3  arcsec,  just  as  we  have  weighted  the  other  astrom¬ 
etry  from  MKO.  These  simulated  observations  reduce  the  major 
axis  of  the  b-plane  ellipse  by  28%,  while  reducing  the  minor  axis 
by  only  6%.  The  impact  threat  derived  with  the  additional  simu¬ 
lated  data  is  given  under  the  label  “Future”  in  Table  13.  The  results 
indicate  that  the  unpublished  astrometry  has  the  power  to  dramat¬ 
ically  reduce  the  impact  risk  from  Apophis,  assuming  that  it  does 
not  significantly  alter  the  new  orbital  solution.  This  goes  to  the 
question  of  orbital  stability,  and  fits  over  different  time  intervals 
indicate  clearly  that  the  debiased  fits  are  very  stable  as  new  data 
are  added,  while  the  raw  fits  show  a  steady  walk  (towards  the 
debiased  solution)  in  the  2029  b-plane  as  new  data  are  added. 


Apophis  2029  Opik  b-plane 
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Fig.  19.  Mapping  of  raw  and  debiased  Apophis  orbital  solutions  to  the  2029  Opik  fa- 
plane,  with  3-sigma  confidence  ellipses.  The  debiased  prediction  (upper  ellipse)  is 
at  the  origin  in  this  depiction  and  the  Earth  is  at  (£,£)  ~  (-9500,-48,000)  km.  The 
Valsecchi  circles,  which  lead  to  a  resonant  return  and  subsequent  impact  are  also 
depicted. 
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Table  13 

Apophis  impacts  in  2036  and  2037. 


Case 

b-Plane  semimajor  axis  (km) 

April  2036 

LOV  sigma 

IP 

April  2037 

LOV  sigma 

IP 

Raw  (JPL#140) 

560 

-2.44 

2.2  x  10~5 

4.11 

8.1  x  10-8 

Debiased 

560 

-3.02 

4.5  x  10“6 

3.53 

8.2  x  10~7 

Future3 

400 

-4.2 

8.8  x  108 

4.9 

3.5  x  10-9 

With  Yarkovsky13 
Debiased 

650 

-2.66 

1.1  x  10~5 

3.10 

2.9  x  10~6 

Future3 

520 

-3.4 

1.7  x  10-6 

3.9 

2.1  x  10-7 

a  Assumes  debiased  solution  is  unaltered  by  addition  of  48  observations  described  in  text. 
b  Assumes  Gaussian  dispersions  of  10  x  10  4AU/Myi  . 


So  far  we  have  not  considered  nongravitational  accelerations  in 
the  trajectory  of  Apophis,  but  we  are  reaching  the  point  where  this 
is  starting  to  become  necessary.  Chesley  (2006)  first  pointed  to  the 
future  importance  of  the  Yarkovsky  effect  (Bottke  et  al.,  2006)  in 
assessing  the  threat  posed  by  Apophis.  Later,  Giorgini  et  al. 
(2008)  conducted  a  thorough  review  of  the  prediction  uncertain¬ 
ties  associated  with  the  Apophis  dynamical  model,  showing  that 
the  Yarkovsky  effect  is  indeed  the  dominant  source  of  dynamical 
modeling  error,  with  potentially  hundreds  of  kilometers  of  offset 
in  the  2029  b-plane  due  to  Yarkovsky  accelerations,  which  is  com¬ 
parable  to  the  current  b-plane  uncertainty.  The  problem  at  hand  is 
that,  lacking  information  on  the  obliquity  of  the  Apophis  spin  axis, 
we  do  not  even  know  whether  the  Yarkovsky  effect  is  causing  the 
semimajor  axis  to  increase  or  decrease.  Chesley  (2006)  found  that 
the  semimajor  axis  drift  rate  da/dt  for  Apophis  should  be  in  the 
range  of  ±15  x  10~4AU/Myr,  but  later  information  suggests  that 
the  object  may  be  somewhat  smaller  than  originally  assumed  (Del- 
bo  et  al.,  2007),  thus  increasing  the  effect.  In  order  to  have  a  rudi¬ 
mentary  understanding  of  how  the  Yarkovsky  effect  might  affect 
the  impact  risk  we  take  a  Gaussian  distribution  on  da/dt  with 
one-sigma  at  10  x  10~4  AU/Myr  and  obtain  a  revised  b-plane  foot¬ 
print.  With  this  dispersion  due  to  Yarkovsky,  the  b-plane  semima¬ 
jor  axis  grows  by  a  modest  amount  as  indicated  in  Table  13,  which 
moves  the  keyholes  closer  to  nominal,  in  a  relative  sense,  thereby 
increasing  the  impact  probability. 

The  Yarkovsky  dispersions  assumed  here  are  unacceptably  sim¬ 
plistic,  and  the  computed  impact  probabilities  suffer  accordingly. 
But  the  exercise  serves  to  emphasize  the  fact  that  Yarkovsky  per¬ 
turbations  have  become  the  key  source  of  uncertainty  in  assessing 
the  hazard  posed  by  Apophis.  We  intend  to  revisit  this  question  in 
the  near  future  with  a  more  rigorous  analysis  that  includes  addi¬ 
tional  astrometry  and  an  empirical  distribution  on  da/dt. 

8.  Conclusions 

We  have  demonstrated  that  our  debiasing  algorithm  signifi¬ 
cantly  reduces  the  biases  in  the  observational  residuals.  We  have 
also  seen  that  only  closely-spaced  observations  of  the  same  object 
made  by  the  same  observatory  remain  correlated.  However,  these 
remaining  correlations  are  significant;  and  an  accurate  observa¬ 
tional  error  model  must  account  for  them. 

The  approach  demonstrated  in  Section  6.1,  while  an  improve¬ 
ment  over  conventional  techniques,  is  sub-optimal,  since  it  essen¬ 
tially  doubles  our  best  estimate  of  the  observational  uncertainties, 
in  order  to  compensate  for  the  unmodeled  correlations. 

Additionally,  we  note  that  combining  observations  from  exist¬ 
ing  surveys  with  the  forthcoming  highly-precise  astrometric  data 
from  PanSTARRS  in  the  same  least-squares  solutions  will  require 
the  use  of  realistic  observational  uncertainties;  simply  assuming 
that  all  observatories  have  the  same  observational  uncertainties 
will  eliminate  any  improvements  that  might  have  resulted  from 
the  more  precise  PanSTARRS  data. 


Our  intent,  therefore,  is  to  use  the  debiased  residual  statistics  to 
develop  observatory-specific,  time-based  models  of  observational 
uncertainties  and  correlations,  as  described  by  Carpino  et  al. 
(2003).  The  resulting  best-fit  orbits  should  allow  us  to  further  re¬ 
fine  our  estimates  of  asteroid  masses,  densities,  and  porosities,  im¬ 
prove  detection  of  Yarkovsky  effects,  and  yield  more  reliable 
collision  analyses. 

The  same  principles  could  also  be  applied  to  develop  a  statisti¬ 
cal  model  of  asteroid  photometry,  potentially  yielding  more  pre¬ 
cise  shape  models  and  rotational  periods.  Combined  with  the 
observational  error  model,  we  should  be  able  to  glean  the  maxi¬ 
mum  information  as  to  asteroid  physical  compositions,  structures, 
and  dynamical  states  from  the  observational  database. 

Note  added  in  proof 

This  analysis  of  the  Apophis  impact  probability  is  included  pri¬ 
marily  as  an  indication  of  the  potential  effect  of  debiased  astrom¬ 
etry.  The  details  of  the  impact  assessment  have  been  updated  with 
new  astrometry  and  a  complete  consideration  of  Yarkovsky  accel¬ 
erations  by  Chesley  et  al.  (2009). 

Acknowledgments 

We  are  grateful  to  Tim  Spahr  and  Gareth  Williams  of  the  Minor 
Planet  Center  at  Harvard  Univ.,  who  provided  the  astrometric 
catalog  information,  which  was  vital  for  this  effort  to  proceed. 
Numerous  observers  and  observing  programs  provided  valuable 
information  beyond  what  is  available  from  the  MPC  data  files.  In 
particular,  we  thank  Bob  McMillan  (Spacewatch),  Ed  Beshore  (Cat¬ 
alina  Sky  Survey),  Scott  Stuart  (LINEAR),  Bruce  Koehn  (LONEOS) 
and  Ken  Lawrence  (NEAT)  for  their  cooperation.  Finally,  we  are  sin¬ 
cerely  grateful  to  Andrea  Milani,  for  his  expertise  and  help  at  each 
step  of  this  effort.  This  publication  makes  use  of  data  products  from 
the  Two  Micron  All  Sky  Survey,  which  is  a  joint  project  of  the  Uni¬ 
versity  of  Massachusetts  and  the  Infrared  Processing  and  Analysis 
Center/California  Institute  of  Technology,  funded  by  the  National 
Aeronautics  and  Space  Administration  and  the  National  Science 
Foundation.  Some  of  the  results  in  this  paper  have  been  derived 
using  the  HEALPix  package  (Gorski  et  al.,  2005).  This  research 
was  conducted  in  part  at  the  Jet  Propulsion  Laboratory,  California 
Institute  of  Technology,  under  a  contract  with  the  National  Aero¬ 
nautics  and  Space  Administration. 

References 

Arias,  E.F.,  Chariot,  P.,  Feissel,  M.,  Lestrade,  J.f  1995.  The  extragalactic  reference 
system  of  the  International  Earth  Rotation  Service,  ICRS.  Astron.  Astrophys.  303, 
604-608. 

Baer,  J.,  Milani,  A.,  Chesley,  S.,  Matson,  R.D.,  2008.  An  observational  error  model,  and 
application  to  asteroid  mass  determination.  Bull.  Am.  Astron.  Soc.  40,  493. 
Bottke  Jr.,  W.F.,  Vokrouhlicky,  D.,  Rubincam,  D.P.,  Nesvorny,  D.,  2006.  The  Yarkovsky 
and  Yorp  effects:  Implications  for  asteroid  dynamics.  Annu.  Rev.  Earth  Planet. 
Sci.  34,  157-191. 


S.R.  Chesley  et  al./ Icarus  210  (2010)  158-181 


181 


Carpino,  M.,  Milani,  A.,  Chesley,  S.R.,  2003.  Error  statistics  of  asteroid  optical 
astrometric  observations.  Icarus  166,  248-270. 

Chesley,  S.R.,  2006.  Potential  impact  detection  for  near-Earth  asteroids:  The  case  of 
99942  Apophis  (2004  MN  4).  In:  Daniela,  L.,  Sylvio  Ferraz,  ML,  Angel,  F.J.  (Eds.), 
Asteroids,  Comets,  Meteors,  IAU  Symposium,  vol.  229,  pp.  215-228. 

Chesley,  S.  et  al.,  2009.  An  Updated  Assessment  of  the  Impact  Threat  from  99942 
Apophis.  Div.  Planetary  Sciences  Mtg.,  October. 

da  Silva  Neto,  D.N.,  Andrei,  A.H.,  Assafin,  M.,  Vieira  Martins,  R.,  2005.  Investigation 
on  the  southern  part  of  the  high  density  astrometric  catalogs  USNO  B1.0, 2MASS 
and  UCAC2.  Astron.  Astrophys.  429,  739-745. 

Delbo,  M.,  Cellino,  A.,  Tedesco,  E.F.,  2007.  Albedo  and  size  determination  of 
potentially  hazardous  asteroids:  (99942)  Apophis.  Icarus  188,  266-269. 

Gilmore,  A.C.,  and  1 1  colleagues,  2004. 2004  MN4.  Minor  Planet  Electronic  Circulars, 
p.  25. 

Giorgini,  J.D.,  Benner,  L.A.M.,  Ostro,  S.J.,  Nolan,  M.C.,  Busch,  M.W.,  2008.  Predicting 
the  Earth  encounters  of  (99942)  Apophis.  Icarus  193,  1-19. 

Gleason,  A.E.,  Larsen,  J.A.,  Descour,  A.S.,  Williams,  G.V.,  2004.  2004  MN4.  Minor 
Planet  Electronic  Circulars,  p.  70. 

Gorski,  K.M.,  Hivon,  E„  Banday,  A.J.,  Wandelt,  B.D.,  Hansen,  F.K.,  Reinecke,  M., 
Bartelmann,  M.,  2005.  HEALPix:  A  framework  for  high-resolution  discretization 
and  fast  analysis  of  data  distributed  on  the  sphere.  Astrophys.  J.  622,  759-771. 

Hog,  E.,  Fabricius,  C.,  Makarov,  V.V.,  Urban,  S.,  Corbin,  T.,  Wycoff,  G.,  Bastian,  U., 
Schwekendiek,  P.,  Wicenec,  A.,  2000.  The  Tycho-2  catalogue  of  the  2.5  million 
brightest  stars.  Astron.  Astrophys.  355,  L27-L30. 

Jedicke,  R.,  and  the  Pan-Starrs  Collaboration,  2004.  The  Pan-STARRS  Solar  System 
survey.  In:  35th  COSPAR  Scientific  Assembly,  COSPAR,  Plenary  Meeting,  vol.  35, 
p.  999. 

Kaasalainen,  M.,  Torppa,  J.,  Piironen,  J.,  2002.  Models  of  twenty  asteroids  from 
photometric  data.  Icarus  159,  369-395. 

Milani,  A.,  Chesley,  S.R.,  Chodas,  P.W.,  Valsecchi,  G.B.,  2002.  Asteroid  close 
approaches:  Analysis  and  potential  impact  detection.  Asteroids  III,  55-69. 

Milani,  A.,  Chesley,  S.R.,  Sansaturio,  M.E.,  Tommei,  G.,  Valsecchi,  G.B.,  2005. 
Nonlinear  impact  monitoring:  Line  of  variation  searches  for  impactors.  Icarus 
173,  362-384. 

Milani,  A.,  Chesley,  S.R.,  Valsecchi,  G.B.,  1999.  Close  approaches  of  Asteroid 
1999ANi0:  Resonant  and  non-resonant  returns.  Astron.  Astrophys.  346,  L65-L68. 


Millis,  R.L.,  and  10  colleagues,  1981.  The  diameter  of  Juno  from  its  occultation  of  AG 
+  0  deg  1022.  Astron.  J.  86,  306-313. 

Monet,  D„  1998.  README  File  for  the  USNO-A2  Catalogs,  <http://www.usno. 
navy.mil/USNO/astrometry/optical-IR-prod/icas/the-pmm/readme.v20>. 

Monet,  D.,  Canzian,  B.,  Harris,  H.,  Reid,  N.,  Rhodes,  A.,  Sell,  S.,  1998.  The  PMM  USNO- 
A1.0  Catalogue  (Monet  1997).  VizieR  Online  Data  Catalog,  p.  1243. 

Monet,  D.G.,  and  28  colleagues,  2003.  The  USNO-B  catalog.  Astron.  J.  125,  984- 
993. 

Pitjeva,  E.V.,  2005.  High-precision  ephemerides  of  planets  -  EPM  and  determination 
of  some  astronomical  constants.  Solar  Syst.  Res.  39,  176-186. 

Skrutskie,  M.F.,  and  30  colleagues,  2006.  The  two  micron  all  sky  survey  (2MASS). 
Astron.  J.  131,  1163-1183. 

Stokes,  G.H.,  Evans,  J.B.,  Viggh,  H.E.M.,  Shelly,  F.C.,  Pearce,  E.C.,  2000.  Lincoln  near- 
Earth  asteroid  program  (LINEAR).  Icarus  148,  21-28. 

Tedesco,  E.F.,  Noah,  P.V.,  Noah,  M.,  Price,  S.D.,  2002.  The  supplemental  IRAS  minor 
planet  survey.  Astron.  J.  123,  1056-1085. 

Tholen,  D.J.,  Bernardi,  F.,  Micheli,  M.,  2008.  Where  is  Apophis?  Bull.  Am.  Astron.  Soc. 
40,  434. 

Thomas,  P.C.,  Belton,  M.J.S.,  Carcich,  B.,  Chapman,  C.R.,  Davies,  M.E.,  Sullivan,  R., 
Veverka,  J.,  1996.  The  shape  of  Ida.  Icarus  120,  20-32. 

Thomas,  P.C.,  Parker,  J.W.,  McFadden,  L.A.,  Russell,  C.T.,  Stern,  S.A.,  Sykes,  M.V., 
Young,  E.F.,  2005.  Differentiation  of  the  Asteroid  Ceres  as  revealed  by  its  shape. 
Nature  437,  224-226. 

Torppa,  J.,  Kaasalainen,  M.,  Michalowski,  T.,  Kwiatkowski,  T.,  Kryszczynska,  A., 
Denchev,  P.,  Kowalski,  R.,  2003.  Shapes  and  rotational  properties  of  thirty 
asteroids  from  photometric  data.  Icarus  164,  346-383. 

Williams,  G.,  2009.  IAU  Minor  Planet  Center,  <http://www.cfa.harvard.edu/iau/ 
mpc.html>. 

Valsecchi,  G.B.,  Milani,  A.,  Gronchi,  G.F.,  Chesley,  S.R.,  2003.  Resonant  returns 
to  close  approaches:  Analytical  theory.  Astron.  Astrophys.  408,  1179- 
1196. 

Zacharias,  N.,  Urban,  S.E.,  Zacharias,  M.I.,  Wycoff,  G.L.,  Hall,  D.M.,  Monet,  D.G., 
Rafferty,  T.J.,  2004.  The  second  US  Naval  Observatory  CCD  astrograph  catalog 
(UCAC2).  Astron.  J.  127,  3043-3059. 

Zacharias,  N.,  and  29  colleagues,  2004.  The  Third  US  Naval  Observatory  CCD 
Astrograph  Catalog  (UCAC3).  ArXiv  e-prints,  1003.2136. 


